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A&stract: 


A comparison of Sky lab and Landsat Cformerly ERTS) images of the 
area around Goldfield, Hev» , shows that for photogeologic interpretation 
of basin-and-range terrain the Skylab images are generally superior. 

The images provide geological data restricted to surficial or 
morphological features such as allwium, playa deposits, bedrock-massif 
borders, lineaments, and, to an extent, alter ed-bedrock areas. Bedrock 
map units or lithologies must be inferred; contacts and actual faults 
cannot be unequivocally recognized. The spatial resolution of the 
Landsat images renders them suitable for small-scale morphological 
studies and lineament analyses. 

The Skylab normal-color stereophotos are an excellent source for 
preparing provisional or preliminary maps of basin-and-range terrains; 
maps derived from these photos would best show surficial units, borders 
of major physiographic features, possible faults and contacts, and very 
broadly inferred lithologic units. These stereophotos have potential 
ll?nited use in the search for alteration zones. Areas of altered rock 
may also be accurately recognized in these photos. 

For many features, spectral-band images showing the greatest 
contrast give the most photo geologic information. Individual-band 
images are not as useful as composite images or normal-color images. 

The value of photogeologic interpretation depends heavily on the 
quality of the images used. For lineament analysis a certain amount 
of graininess or misregistration seems to improve recognition, however. 
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T6.e linages used by us in this study were those couHsercially 
available from the EROS Data Center, Sioux Fails, S. D. Of the images 
available to us for this test area, the most useful for photogeologic 
analysis in the Skylah system are the normal-color stereopair and the 
single-band 0.5-0^ 6 pm image. The most useful images in the Landsat 
system are the false-color IR image and individual bands 5 and 7. The 
Skylab normal-color stereopair is especially useful for recognition 
and mapping of altered areas. When compared with procedures involving 
image enhancement, analyses based on relationships of pattern, color, 
topography, and alluvium provide quick and perhaps surprisingly 
accurate identification of altered areas. 
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Introduction 


This study was made to determine in what ways haiidsat Cfoinnerly 
ERTS) images differ from Skylab images as tools for photogeologic 
interpretation of the area around Goldfield, Nev. Our object was to 
decide which of the two imaging systems provides more geological 
inforination on geologic units, structural features, lineament patterns, 
and landforms. 

As Rowan and others (1974) have amply demonstrated, interpretation 
of geologic features in this region from Landsat images is potentially 
of great utility and precision* Accurate photogeologic interpretation 
of the whole imaged area or of selected features will save considerable 
time and expense where field work is especially arduous and where 
access is limited. The main value of unenhanced high-altitude images, 
such as Skylab photos, is in reconnaissance; they provide an acctarate 
basis on which the geologic relationships of major physiographic 
features may be reasonably inferred, so that groundirork may be 
efficiently directed to areas of particular interest* Certain problems 
must be coxintenanced, however. In this study these included degraded 
resolution and poor color balance, probably due to variations introduced 
by the fact that the prints and negatives were several generations 
removed from the original images. Lack of good ground control and 
suitable geologic data for comparison with the interpretations are 
problems that also influenced our conclusions. 

Bsspite these restraints, a precise photogeologic map of the study 
area can be produced from these images, although it must be essentially 
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a pTiysiograplixc-surficial map. The photogeologic-interpretation 
product is rated here in relation to spectral and spatial resolution 
of the images and is compared with published geologic maps to arrive at 
a general evaluation of accuracy, scope of application, and ease of 
recognition of pertinent geologic features. 

Study area, geologic setting, and images 

The study area, shown in landsat frame E 1324-18005, dated June 
1973, includes parts of southwest Nevada and southeast California 
Cfig. i) . The portion of California (Inyo County) includes parts of 
the Inyo fountains, the Panamint Range, Grapevine Mountains, and Beath 
Valley C^ig. 2); the portion of Nevada (Esmeralda and Nye Counties) 
includes Mud Lake, the Montezuma Range, the Cactus Range, part of the 
Kawich Range, and the. northern part of the Grapevine Mountains (fig. 2). 
The principal town in the area is Goldfield, Nev. This region has 
excellent bedrock exposure and generally clear atmospheric conditions, 
which prevailed at the time the data were recorded. Vegetal cover is 
sparse or absent. Control information for the whole area is available 
in geologic maps of both Nevada and California (Jennings, 1958; 

Strand, 1967; Stewart and Carlson, 1974).. 

Physiographically, the study area is characterized by rugged, 
relatively deeply dissected, fault-block mountain groups. These massifs 
are more-or-less surrounded and partly buried by alluvial fans and 
playas. Elevations range from 11,107 feet (3,387 m) on Mt. Inyo, Calif,, 
to -200 feet (-61 m) in Death Valley, Calif. In most cases drainage 
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is closed or ephetaeral. The physiography is shown in detail on the 
following 2*^ quadrangles, at scale of 1:250,000; Death Valley, Calif., 
Nev,; Goldfield, Nev., Calif.; and Tonopah, Nev. 

The general geology is that of the Basin-and-Range Province: 
chiefly clastic an ^ carbonate units variously faulted and tilted and 
locally intruded or overlain by volcanic rocks. Rocks are present to 
lesser extent, especially near the Calif ornia-Nevada border. The 
massifs are mantled by coalescing alluvial fans and talus cones. These 
feed down to braided stream washes, which lead to playas, ephemeral 
lake bottoms floored with clay and evaporites. 

The images used from Skylab are the Sl90~A 70-mm photographs in 
the following spectral bandwidths: 0.5-0. 6, 0.6-0. 7, 0.7-Q.S, and 

0. 8-0.9 pm, and a color-IR print. A normal-color stereopair, in which 
the Skylab-Landsat image overlap delimits the study area, was also 
Included. The corresponding Landsat images are band 4 CO. 5-0, 6 pm), 
band 5 CO* 6-0. 7 pm), band 6 CO.7-0.8 pm), and band 7 CO, 8-1.1 pm), and 
a commercially available color-IR-image composite of bands 4, 5, and 
7. Scales of the images are between 1:500,000 and 1:1,000,000. All 
of the images are paper prints. 

Method 

Our method of comparing the Skylab and the Landsat images for 
this study area is basically an item-by-item comparison keyed to data 
presented in the geologic maps of the area. The items of geological 
interest are those that can be directly interpreted from the unenhanced 
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fmags, defined by a boundary line, and drawn on a map at fcbe scale of 
tbe image. Such items are functions of morphology, structure, and 
material composition. Because only morphology is directly observable — 
structure and composition can only be iriferred from the images— we 
ought logically to choose landform patterns as the geologic units 
rather than lithologies or even formations. This approach is 
especially practical where a good stereoscopic model is available. 
Otherwise, we are restricted to one- and two-dimensional elements of 
landforms; i.e. , linear elements, represented hy such features as 
gorge bottoms, stream segments, scarp edges, ridge crests and other 
slope divides, and straight boundaries between surfaces of different 
texture, and planar elements of distinct textural and tonal character, 
such as various kinds of slopes and shapes. Morphological units 
Cfeatures made of two or more slopes) are of two general kinds: 
erosional landforms (iuselbergs and massifs) and depositional landfoms 
Cvolcanic flows and cones, playas, and alluvial fans). The kinds of 
photogeologic features considered in this study are listed in Table 2, 
along with appropriate features shorn in geologic maps. 

The features recognized and defined in the images were drawn on 
an overlay map. Because of uncertainties in the mapped geology, 
devising an objective statist‘s cal test to evaluate these data, except 
for the lineament data, was impossible. The lineament analysis was 
done by one operator. To minimize increased recognition of features 
due to familiarity with the image pattern and, hence, uncontrolled 
bias, lineaments were recorded over a three-week period and each 
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image was accorded the same amount of time during analysis. 

STevertheless, a "learning" effect developed and a rough, correlation ; 

• j 

exists (Table 1) between the order of mapping and the number of i 

i 

lineaments per image correlated with mapped faults compiled from the i 

] 

j 

State maps. . ] 

j 

I 

Photogeologic observations j 

The following paragraphs summarize comparisons of each kind of 
mapped feature observed in the images. An evaluation of the images 
with respect to each feature is given. i 

Playas- Recognition and shape discrimination of playas VTas better i 

in the normal-color Skylab images than in the Landsat color composite, 
because of the better spatial resolution, the color contrast, and the j 

st er eomo deling . 

Twenty playas were recognized in each of the Skylab images (fig. 3). | 

of the single-band photos^ the 0,5-0. 6 band image provides the best | 

discrimination because of high tonal contrast in this bands recognition i 

in other photos was difficult because of lower contrast and excessive | 

graininess. Thirty-two playas were Identified in each Landsat image I' 

(fig. 4); of these images, bands 4 and 5 provided the greatest ease of | 

% 

recognition because of the greater tonal contrast. On the other hand, 
tonal variations in the larger playa surfaces are shown in bands 6 | 

and 7. These variations indicate subtle contrasts in texture and, 
probably, moisture content. Ail 20 playas recognized in the Skylab 
images were among those recognized in the Landsat frames. 
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i'ifteen. playas are shown on geologic maps of the. study area. It 
was detertEixned that 27 of the 32 Landsat-identified playas were 
spurious: some highly reflecting material other than playas led to 

misidentification. Of the 20 Sky lab-identified playas, 6 were 
misidentified for the same reason. One playa shown on the geologic 
map was not recognised in these photos. Clearly, the Sky lab photos 
are a superior source for identifying playas in basin-and-range terrain. 

Alluviunt - Source maps for the California State map distinguish 
among different kinds of alluvium, but the compiled maps show it as a 
single general unit. The images all allow recognition of stream-laid 
alluvium versus cone, talus, or fan deposits Cfigs. 3 and 4), and the 
relationships correspond with those mapped on the California 
quadrangles. In general, the Sky lab photos provide the better 
discrimination because of superior spatial resolution. Spectrally, 
the two systems seem to respond differently: the Landsat hand-5 image 

is the best of the Landsat bands for distinguishing alluvium, hut the 
Sky lab frame 0. 5-0.6 is the best of the Sky lab hands. The quality of 
the 0, 6-0.7 fra.me in the Sky lab series (which corresponds to Landsat 
band 5) is poor, probably because of a processing irregularity; however, 
the cause is uncertain as a valid hand-to-hand comparison between the 
two systems was not possible. Of all the images, the Skylah 
normal-color stereopair is the best, because of excellent chromatic 
contrast and the stereomodel. 

Bedrock areas - Bedrock boundaries drawn from the Landsat images 
(fig. 4) are more deeply embayed and circumscribe smaller areas than 


8 


do those drawn from the Sky lab photos Cfig* 3), owing to the low tonal 
contrast between much of the alluvium Ctalus deposits) and the adjacent 
bedrock. Because much of the distinction depends on texture, the 
relatively low spatial resolution of the Landsat images puts them at 
a disadvantage. Thus, a conservative estimate of bedrock areas tends 
to err in favor of alluvium. The boundaries observed from the Sfcylab 
color photos show the best correspondence tri.th mapped contacts. The 
borders were recognized in these photos by the additional tonal and 
color contrast and by the stereorelief. The stereorelief is especially 
helpful bt^'^ause of the sharp morphological distinction beiiween bedrock 
and alluvium. However, bedrock with very shallow fan cover (pediments), 
bedrock with albedo similar to adjacent fan deposits, rotten or 
disaggregated bedrock with surface texture similar to alluvium, and 
small, scattered outcrops below the spatial resolution of the sensor 
were included as alluvium. In the California portion of the area, for 
example, younger nonmarine units are not distinguishable from older 
nonmarine units, because these units have spectral characteristics 
similar to those of recent alluvium and have no significant physiographic 
expression (Jennings, 1958). 

Of the Landsat images, the false-color composite and band 5 offer 
the best tonal contrast for determining outcrop boundaries. Band 5 
has the better tonal contrast, but it suffers from striping, Ey far 
the best Images for bedrock definition are the Skylah normal-color 
stereopair. The color-IR print offers very poor tonal contrast, 
and the individual band images are too grainy or have low tonal contrast. 
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Bedrock, units - In the Nevada portion of the study area, geologic 
units Interpreted from Both Skylah and Landsat images show no 
relationship to units presented on the State map, because the map units 
do not represent unique, spectrally distinctive lithologies. In the 
California portion, however, the Skylab photos provide surprisingly 
good identification of contacts but not lithologies. Photogeologic 
units are distinguished according to texture, tone, and color. The 
volcanics are particularly easy to recognise in both Landsat and 
Skylab photos, because of their distinctive morphology and dark tone 
and color (figs. 3 and 5). Band 4, Landsat, is an inferior image for 
this purpose, however; band 0,5“0,6 is the best of the Skylab group. 

Not all recent volcanics shown on the maps were recognized in the 
images because of variations in image quality. Some features 
interpreted as volcanics cannot be confirmed, because they fall in 
unmapped areas on the geologic maps; others appear to have been 
incorrectly identified, probably because they are located near the 
edges of images, where resolution fall-off occurs. 

Contacts in the Grapevine Mountains are particularly well 
distinguished in the Skylab normal-color images (fig. 3) . In the area 
south of Dry Mountain, at the west end of the Cottonwood Mountains, 
Precam.br ian and Pennsylvanian units are well defined. In the Inyo 
Mountains, along the west side of Eureka Valley, two contacts between 
Cambrian and granite units are identified. And at the north end of 
the Last Chance Range, the Furnace Creek fault zone and a contact 
between granite and Cambrian rock is well defined. The recognized 
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contacts separata map units of broadly different lithologies- and 
textural and (based on the stereomodel) -morphological expressions. 

Low tonal contrast hampered attempts to identify geologic imits 
from the Landsat images. Volcanics are fairly well deflr4.ed on the 
images Cfig- 5), but the shapes and locations of many volcanics are 
not accurately determined. 

“For these data, the Sky lab normal-color photos have the best 
contrast, and the color-IR photo, the worst. Spectral resolution in 
the individual band images is too low to be useful for discriminating 
bedrock units both in the Sky lab and in the Landsat images. 

Lineaments - Most of the lineaments interpreted from the Skylah 
Landsat images do not correspond to mapped faults, chiefly because 
the mapped faults are short segments within massifs- or along borders 
of massifs and alluvium, whereas the lineaments include features 
outside of massifs. Lineaments that coincided with mapped faults were 
counted in each of the 11 images. Results are shown in Table 1. 

The Sky lab black-and-white images appear to be more useful than 
the Landsat images for mapping lineaments, probably because of the 
higher contrast, which produces sharp grain edges on the Sky lab images, 
and because of the higher resolution of the Sky lab images. However, 
the Landsat color-IR image appears to be superior for this purpose to 
either the Sky lab color-IR or normal-color photos, possibly because 
of the slight edge enhancement caused by imperfect registration of the 
three bands used to produce the Landsat false— color composite. 

General summary - The Sky lab stereophotos are recommended for 
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use in Tnaklng surficial geologic maps and for making preliminary 
geologic maps in arid terrains, in which the map units include areas 
of bedrock exposure, playas, alluvial fans, stream-laid alluvium, and 
lineament patterns. In other areas, dune deposits, swamp deposits, and 
glacial deposits may he defined. Broadly defined lithologic contacts 
may be outlined from these photos, and, to an extent, preliminary 
identification of altered areas in arid regions may be achieved. 

Interpretations and restrictions 

Our interpretations depended directly on the spatial and spectral 
characteristics of the images used here. The Skylab photos have a 
relatively high spatial resolution C30 to 75 m), whereas the Landsat 
resolution is lower (80 m) . Because of the better spatial resolution 
and the excellent stereomodel offered by the Skylab photos, they are 
superior sources of physiographic data. The high resolution is not 
necessarily of value for recognizing the extent and position of 
lineaments and other features that compose a large number of resolution 
elements, neither image system provides a superiority in shadow 
enhancement, as there is only a 4° difference in snn-elevation angles 
between the images (62° for Landsat, 66° for Skylab) . Likewise, where 
similar handwidths are covered, spectral resolution is similar in the 
multiband images for both systems. The Skylab normal-color photos 
offer the advantage of good resolution combined viith the erlianced 
visual discrimination of color tones relative to gray tones. The 
spectral contrast of the Skylab color -IR photo, on the other hand, is 
lowj and the photo is not particularly useful. 
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Despite the general utility of Skylab color photos, both Shylah 
and Landsat data have serious limitations. Our interpretations of the 
images, or of the photo geologic units, led to a geomorphic synthesis 
of the study area. We found that, though useful, the geologic 
information contained in each kind of image was not suitable for 
making a geologic map of the area and was only marginally useful for 
improving the information already contained in the maps for two 
reasons: Ci) conventional geologic maps- depict lithologic units, which 

are not unequivocally identifiable in Skylab or Landsat images, and 
(2) most lineaments, to be documented as fault traces, must be 
field-checked. Alteration effects are best displayed in a 
special-purpose map ordinarily made in the interests of economic geology. 
The Nevada State map, for example, is a preliminary compilation; the 
map units are essentially lithologic units grouped according to broad 
compositional and textural distinctions; e.g,, "Intrusive rocks: 
aphanitic, porphyritic, and granitic rocks ranging in composition from 
diorite to granite," or "Ash-flow tuffs, rhyolite flows, and shallow 
ir.trusive rocks." The correlation of units based on wide ranges of 
origin, grain size, and bulk composition with units based on spectral 
reflectance is virtually impossible. 

An attempt was made to evaluate the geologic data control in the 
regional maps. An examination of the Mina quadrangle ^Ferguson and 
others, 1954), a good example of detailed mapping in the study area, 
revealed the same inadequacies. A sparsity of attitude symbols in 
the central and southwest quarter of the map indicates poor control 
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for structural interpretation j fault napping is inadequate tftroughout 
nost of tlie map area. Map units embrace abrupt facies changes and 
changes of metamorphic grade, such as "sandstone (tuffaceous in part), 
shale, calcareous marl, limestone, diatomite and local basal 
conglomerate. Locally interbedded rhyolite tuff is present." 

Crystalline units, based on a broad range of compositions, are 
lumped. Intrusive rock units, for example, contain rocks ranging 
from granite to gabbro to serpentine. 

Tvro source maps in the California section were examined (Blackcap 
Ifountain quadrangle, Bateman, 1965; Blanco Mountain quadrangle, Nelson,. 
1966) , These maps show units grouped according to origin and age 
(i.e., marine vs nonmarine) but not according to lithology. In some 
areas, units mapped by reconnaissance are sharply juxtaposed with 
units mapped by detailed work, with blank "unmapped" borders. Both of 
these maps show units and structure in excellent detail. The problem 
here, however, is too fine a division of units. Bor example, in the 
Blackcap Mountain quadrangle, there are four units of quartz monzonite 
and five units of grandodiorits. Individual map units of this sort 
were indistinguishable in the Skylab and. Land sat images; although 
there is a broad matching of contacts between widely different 
lithologies, it is impossible to distinguish one granodiorite or one 
quartz monzonite from another in the images. 

Rowan and others (1974) , faced with the same problem but with 
considerably superior image data, compared an "excellent color 
photograph" obtained from the S-190 photographic experiment with a 
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map showing the distrihutiou of the main rock types. They concluded 
that the appearance of the rocks in the photo is dominated hy albedo, 
^ich is not a reliable guide to rock type; the characteristic colors 
of the 7:ock units are too muted to be consistently shown on small-scale 
orbital photographs. Further, except for a few small limonitic 
altered areas southeast of Stonewall Mountain, Nevada, the mineralized 
areas are not distinctive in the Skylab photograph. Rowan and others 
noted that rock- type discriminations on the standard Landsat MS'S images 
are severely limited, even if only a two-component classification system 
of mafic and felsic rocks is used; spectral-reflectance differences 
among rock types and between altered and unaltered rocks are generally 
too small to be detected by visual comparison of the MSS images or 
through analysis of color-IR composites. 

The comparative study of Skylab and Landsat images for the purpose 
of lineament analysis does not make optimum use of the better 
resolution of the Skylab system. If all the resolvable lineaments 
were portrayed, the resulting illustration would show a cloud of lines 
representing lineaments tens to hundreds of metres in length. In 
order to deal with a manageable amount of data, lineaments in the 
analysis were limited to an arbitrary minimum length of 3 km, 

A systematic study of lineaments shows that the percentage of 
lineaments that coincide with mapped faults differs according to the 
images used for data. This variation is due chiefly to the effects 
of contrast, photographic grain, and edge enhancement. Visual 
Inspection of the photos shows that the Skylab photos are generally 
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Higher in contrast than the Landaat images. So many steps of photo 
processing have occurred between the original data and the images used 
for the analysis, that it is impossible to determine whether the higher 
contrast is due to inherent properties of the systems or is a product 
of the photographic-reproduction process. However, the Landsat prints 
have to be considered as an example of what the user can expect in the 
way of working materials. A subjective judgment is that, in general, 
the Skylah prints used in this study appear to he grainier than the 
Landsat images. Sharp grain boundaries in the Skylab photos may 
actually enhance the contrast associated with most lineaments. Grain 
boundaries in the Landsat images appear to be soft. Edge enhancement 
is important in the Landsat false-color composite, where a slight 
misregistration of three separate images is present. A close examination 
of this print shows that the misregistration is due to both a shift in 
registration of the bands and a slight scale change from band to band. 
This effect of broadening individual lines seems to facilitate 
recognition of a greater number of lineaments than can he obseived on 
the individual band images. 

Although spectral resolution figures in the recognition of 
lineaments, this phenomenon is also difficult to evaluate. Because of 
decreased scattering at longer wavelengths and hence better atmospheric 
penetration and concomitant blacker shadows, we would expect more 
lineaments to be recognized on longer wavelength images. In fact, 
this may be reflected in the Landsat data in Table 1, although the 
small increase in lineament recognition with wavelength may not be 
statistically significant. 
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FurtTiisr potentials for pliotogaologic analysis 

Practical applications of satellite-lasaage data have, yet to be 
fully explored. Rowan and others C1974) have made impressive advances 
using Landsat MSS data in the search for mineralized areas in Nevada, 
but with respect to using unenhanced images and photogeologic 
interpretations for this purpose, they are not enthusiastic: 
"...preliminary evaluation suggests that small-scale photographs such 
as the Skylab example, although useful for morphological and structural 
studies, are not adequate for detecting and mapping mineralized areas" 
Cp. 12), Their objections, however, are based on the use of color as 
the sole criterion for recognition of alteration; the Skylab stereopboto 
pair available to us provides us with a variety of diagnostic features. 

Critical features in altered areas are related to topography and 
microrelief. McKinstry C1948, p. 223) pointed out a basic relationship 
useful in the search for ore deposits: 

"In the arid basin-and-range country of the 
southwestern United States, pediments are specially 
favorable situations for vein deposits. 1'Jhereas 
the mountains represent in general the more massive 
and resistant rocks, the pediments develop on 
those parts of the rock mass which faulting, 
fracturing, and chemical decay have rendered 
vulnerable to erosion. It is in just these places 
that hydrothermal ore deposits are to be expected, 
not only because zones of structural disturbance 
are hospitable, but because hydrothermal solutions 
soften the rocks and because pyrite, by furnishing 
sulfuric acid, promotes weathering." 

iinothp.-r topographic association occurs where a resistant, unaltered 

rock overlieii altered rock, such as that described by Anderson and 

others (1965) . Based on field work, they suggested a search for 
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concealed mineral deposits in the area between Stonewall Mountain and 
Mountain Helen, and north of Mountain Helen where isolated patches of 
altered rocks are surrounded by the unaltered Thirsty Canyon Tuff. 

This area is characterised by ridges of gray to red-gray dense, 
silicified rocks that are bounded by valleys formed in light-gray, 
pink, and tan, intensely argillized rock. In many places the altered 
areas abut the tuff and probably extend beneath it, 

• "S' 

A search for such altered areas was made on the Skylab normal-color 
stereopair, which partly overlaps the study area of Rowan and others 
(1974). Approximately 50 percent of the altered areas that they 
detected in this restricted area were recognized in the Skylab photos 
(figs. 6a, b). Other areas included rocks similar in color and 
weathering to altered materials. Two sites were spotted at the north 
end of Stonewall Mountain that were not detected on the image of 
Rowan and others; these sites are marked by mines or prospects. Areas 
missed were principally those where cloud cover made interpretation 
very doubtful. Probably more altered areas would have been detected 
had the scene been cloudless. These areas were identified entirely 
on the basis of the criteria listed below; during the search no 
reference was made to the results of Rowan and others. The results 
obtained in this study indicate that, lacking sophisticated 
image-enhancement facilities, areas can be located on the photos that 
are worth investigating on the ground. This approach is quick and 
rough, but it is cost effective and it extends the usefulness of the 
raw Skylab images to more geologists. The use of cloud-free photos 
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printed to maximim sharpness and correct color balance would help 
greatly in this pursuit. 

The criteria used for interpreting altered areas are as follows: 

1) Pattern: areas show a mottled, "stained,” irregular pattern 

of generally light tones. Borders are typically 
irregular but are sharp. 

2) Color; colors are distinct from surrounding rock: typically 

off-white Cclay and silica) , pale tan Climonite, 
sulfate staining) and/or rusty red Chematite) . These 
colors are commonly mottled together. 

3) Topography: areas are adjacent to harder rock; hence they are 

recognizably bedrock areas but are rather low. 

Because of variable resistance of altered products 


4) Alluvium: 


and association with faulting, the areas are | 

topographically rough: highly dissected with a i 

1 

i 

subdued, bumpy, or somewhat deflated appearance. { 

I 

Adjacent hard rock may form an obvious cap over or ^ 

around such areas, or be present as knoblike outliers, | 

1 

As HcKinstry (1948) pointed out, such areas are ^ 

typically adjacent to pediment borders; hence a good ' 

deal of distinctive slope wash is present, streaked ; 

i 

with the typical alteration colors, | 


Conclusions 

Sky lab photographs provide more and better detailed information 
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at scales useful for topical studies and initial detailed mapping 
projects in terrain of the basin-and-range type; Lar.dsat images- provide 
a useful small-scale data base to prepare regional geomorphic-tectonic 
compilations or syntheses -where much geologic mapping has already been 
done. The images are marginally useful in delineating broad areas of 
distinct lithologies. In most cases it is difficult to reconcile 
photogeologic information with mapped geology, because individual 
concepts have gone into making a map and because the definition of a 
formation is not always compatible with a spectrally distinct lithology 
as shown on a remote-sensor image. The Sky lab photographs are 
especially valuable because of the excellent stereomodel they provide. 
On the whole, the most useful photos in the Sky lab system are the 
normal-color stereopair and the single band 0.5-0. 6 pm. 

The Landsat false-color image and individual bands 5 and 7 are the 
most useful in the landsat system for small-scale regional-landform 
and tectonic-pattern analyses. 

Areas of altered rock can be recognized in the Skylab 
normal-color stereophotos using the criteria outlined in the previous 
section. The image-enhancement procedure developed by Rowan and others 
C1974) provides a more abjective, machine-controlled process in which 
spectral information is gained from the images. On the other hand, 
this machine process causes information on topography and texture to 
be lost. This lost Information is essential to the identification of 
altered areas through the use of Skylab photos. To make the most of 
all the information potentially available to the analyst, both 
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procedures should he used in conjunction. Otherwise, the Shylah photos 
alone may provide relatively quick and cheap reconnaissance of altered 
areas . 


21 


References cited 


Anderson, R. E. , Ekren, E. B. , and Healey, D, L, , 1965, EosslBle buried 
mineralized areas in Hye and Esmeralda Counties, Nevada, in 
Geological Survey research 1965: U S. Geol. Survey Prof. Paper- 
525-D, p. D144-150. 

Bateman, P. C. , 1965, Geologic map of the Blackcap Mountain quadrangle, 
Fresno County, California: II. S. Geol. Survey Geol. Quad. Map 

GQ-428, 1:62,500. 

Ferguson, H. G., Muller, S. W. , and Cathcart, S. H. , 1954, Geologic 

map of the Mina quadrangle, Nevada: tf.S. Geol. Survey Geol. Quad. 

Map GQ“45, 1:125,000. 

Jennings, C. W. , 1958, Geologic map of California — Death Valley Sheet: 
California State Dept. Nat. Resources, 1:250,000. 

McKinstry, H. E. , 1948, Mining geology: New York, Prentice-Hall, Inc. , 

680 p. 

Nelson, C. A., 1966, Geologic map of the Blanco Mountain quadrangle, 

Inyo and Mono Counties, California: II. S. Geol. Survey Geol. Quad. 

Map GQ-529, 1:62,500. 

Rowan, L, C., Wetlaufer, P. H. , Goetz, A. F. H. , Billingsly, F. G., 
and Ste-wart, J. H. , 1974, Discrimination of rock types and 
detection of hydro thermally altered areas in south-central 
Nevada by the use of computer-enhanced ERTS images: U.S, Geol. 

Survey Prof, Paper 883, 35 p. 


22 




Table. 1, percent correlation of lineaments and mapped faults in each. 
Sk.ylab and Landsat image, Goldfield area, Nevada 


Skylab S190A images 

Image 

% correlation 

Order of lineament 
mapping 

Number of 

lineaments observed 

0.5 - 0.6 pm 

25.0 

8 

124 

0.6 - 0.7 pm 

21.8 

9 

119 

0.7 - 0.8 pm 

34.0 

10 

94 

0,8 - 0.9 pm 

24.1 

11 

IQB 

normal color 

13.7 

1 

183 

color IR 

17,4 

2 

178 


Landsat iicages 


Image 


% correlation 

Order of lineament 
mapping 

Number of 

lineaments observed 

0.5 - 

0.6 

pm 

C 4 ) 

14.2 

4 

127 

0.6 - 

0.7 

pm 

C 5 ) 

14.8 

5 

168 

0.7 - 

0,8 

Pm 

C 6 ) 

15.4 

6 

175 

0.8 - 

1.1 

pm 

( 7 ) 

17.6 

7 

199 

false- 

■color IR 

22.5 

3 

193 
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Table 2, Criteria for 


Photogeologic 

feature 


1 . lineament 


2. erosional 
landforms 
(bedrock 
areas - 
inselbergs ) 


3, depositiotial 
landforms 
(alluvium; 
stream deposits, 
cone deposits, 
playas) 

4. lithologic units 


recognizing photogeologic 


Criteria 


location and extent of 
aligned tonal or 
topographic features 
(i.e. , streams, gorges, 
alluvium-h edro ck 
contacts, ridges, 
textural boundaries, 
etc.) along straight 
or slightly curved lines 

extent and outline of 
inselbergs and massifs 
based on relief (stereo- 
model) texture, tonal 
contrast, surface (e.g,, 
drainage) patterns, 
shadowing , and color 

extent and outline based 
on texture, tonal and 
color contrast, surface 
patterns, and boundary 
relationships 


delineation and 
interpretation based 
on color, tone, 
texture, boundaries 


features 


Map unit 


mapped and projected 
faults and major 
unit contacts 


all bedrock units 


alluvium, playa 
deposits 


geologic units 
(formations) 


25 



I ^ ^^rsoKi C.ii'Y 



o Jo ^ too //6* 

hito^e’f‘r&S 


F~i^Uf€. ! , Loc&J’ioiri ^ka.p€. o-f 4-hc S'ku^/y 






Genera! phy -f€a''f uf'€^ <5/ //>e 

S’fuc!^ a^rea^ . 





EXPLANATION 



Stream- \^id alluvium 

Fan deposits 

■ Q ^ I Volcanic rocUs-- 

Mostly mafic 


clouds 


Massif areas and 
inselbergs-Lines 
represent distinct 
lithologic contacts 

-- Mapped faults 
coincident with 

lineaments 


•\ ‘ 


Lineament; 
dashed 
where ob 
scure 


STATE 

BNDY. 


scale approx. 


} 

i 


Figures. PhoTOGEOLOGIC map MADE FP.OM SKYLAB NORMAL 

COLOR. 5TEREO IMAGES reproducibility of th£ 

^ IT , ORIGINAL PAGB IS POOR 




EXPLANATION 


Massif areas and insd 


Volcanic rocks-- Mos+I^ mafic 


Stream- laid alluvium 


Lineaments-- 
D^ished where obscure 


STATE 

BNDY 


25 Krn 


SCALE APPROX 


Figure 4. Photogeologic map made from LAMDSAT 

COLOR IMAGE 


K ^ 



Figure 5. Areas of volcamic rocks, inferred to be 

ALL OR MOSTLY MAFIC, MAPPED FROM LAMDSAT 
BAND 7 IMAGE 


reproducibility of the 

ORIGINAL PAGK IS POCK 






Fi<gure 6a. Areas of hydrothermally altered rock. 

MAPPED FROM SkYLAB NORMAL COLOR 
STEREO IMAGES 
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Figure 6b. Areas of altered rock identified by Rowan 

ET.AL.{/?74), COMPARED WITH THOSE IDENTIFIED 
FROM SKYLAB normal color STEREO IMAGES 
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Figure I. IxKation maps of test sites; A. four F.RIM flight 
lines in southern California and Nevada; R. Skylab II and III 
tracks over southwestern California and eastern Nevada. 


ABSTRACT 

High-reflectivity anomalies in the 
near infrared (1.0 to 2.6 nm) have been 
observed on scanner images obtained by 
the Environmental Research Institute 
of Michigan (ERIM) in flints over 
several geologic units in southern Cali- 
fornia and southwestern Nevada and by 
the Skylab 11 mission in flights over 
western Nevada. These anomalies almost 
always occur in rocks that are at least as 
mafic as andesite. The most reasonable 
explanation for the .anomalies is that the 
anomalous materials have a 5-fim or 
finer coating of hematite that was formed 
during cooling and posteruptive fuma- 
rolic steaming. 

INTRODUCTION 

Prior to the launching of Skylab, 
the Environmental Research Institute of 
Michigan (ERIM), in several aircraft 
missions flown for the U.S. Geological 
Survey, used a 12-channel infrared 
scanner to simulate the various channels 
of the Skylab S-192 experiment. (The 
S-192 experiment is an optical mechani- 
cal scanner that operates in 13 bands 
located in visible and infrared regions 
with high atmospheric transmission.) 

Data from both missions showed near 
infrared reflectance anomalies of ande- 
site and basalt in southern California 
and Nevada. 

In September 1972, the ERIM 
12-channel scanner was flown over six 
areas believed to be representative of 
the major geologic units in the proposed 
southern California-southwestern Nevada 
Skylab test site. Although operational 
constraints prohibited surveillance of 
this particular area by Skylab, u geo- 
logically similar test site in western 
Nevada was overflown by the spacecraft 
and similar data on reflectance anomalies 
were obtained. 

The ERIM aircraft scanner was 
flown in a mode that permitted data 
collection in nine channels in the ultra- 
violet, visible, and reflective infrared, 
and three channels in the thermal 
infrared. Only the eight channels in 
the visible and reflective infrared (listed 
on the following page) are discu.s.scd in 
this report. 





Figure 2. ERIM scanner images of the cinder cones northeast of Old Dad Mountain; 

A, channel 6 (0.67 to 0.94 /Jm); B, channel 8 (I.O to 1.4 jL(ni);C, channel 9 (2.0 to 2.6 /Jm). 


ERIM Channel Wavelength (ium) 

1 0.41 to 0.48 

2 0.46 to 0.49 

3 0.50 to 0.54 

4 0.55 to 0.60 

5 0.62 to 0.70 

6 0.67 to 0.94 

8 1.0 to 1.4 

9 2.0 to 2.6 

The Skylab scunner was flown in a 
mode that permitted data collection in 


13 channels. Only two have been used 
in this report. 

Skylab Channel Wavelength (jam) 

7 0.78 to 0.88 

11 1.55 to 1.75 

Both the ERIM te.st sites and the 
Skylab Nevada test site contain most of 
the igneous and sedimentary rock types 
with the exception of ultramafic. Granite, 
andesite, basalt, and limestone are par- 


ticularly common. A location map of 
four of the six ERIM sites, as well as the 
Skylab Nevada test site, are shown in 
Figure 1. 

DATA 

In the images obtained from the 
ERIM missions, many small areas of 
higher or lower reflectivity were noted 
from channel to channel; however, the 
most striking contrasts v./ere observed 
among channels 6, 8, and 9. A typical 
example is illustrated in Figure 2, which 
shows an area just northwest of Old Dad 
Mountain (lat 25° 1 2' N., long 1 1 5°47' 
W.). The area is mapped as Quaternary 
cinder cones underlain by Quaternary 
basalt (Barca, 1966). The 0.67 to 0.94 
Hm channel image is similar to the 
visible channel images; that is, the low 
reflectivity of the cones and subjacent 
basalt is relatively constant below 1 fim. 
The summits of the cinder cones, how- 
ever, exhibit moderately high reflec- 
tivity in channel 8 ( 1 .0 to 1.4 fim) and 
very high reflectivif in channel 9 
(2.0 to 2.6 iim). This type of anomaly 
was also observed in an area east of 
Pisgah Crater and southwest of the 
town of Ludlow (lat 34°40' N., long 
11 6° 14' W.)and in the Cady Mountains 
(lat 34° 52' N., long 1 16°21' W.), as well 
as in a number of places shown in the 
Skylab images of the Nevada test site. 

In the Skylab S-192 images, anomalies 
appear in Skylab channel 11(1 .55 to 
1.75 /am), but not in images at wave- 
lengths shorter than 1 .0 (channel 7, 
0.78 to 0.88 /am. Fig. 3). In both in- 
stances, the anomalies arc almost always 
in volcanic rocks that arc reddish, 
mapped as Tertiary or younger, and at 
least as mafic as andesite (Ross, 1961 ; 
Dibblee, 1966, 1967a, 1967b; Dibblec 
and Bassett, 1966; Strand, 1967; Corn- 
wall, 1972; .Albers and Stewart. 1972); 
however, most of the andesitic and 
more mafic rocks in these areas do not 
show anomalous infrared reflections. 

DATA REDUCTION 

In an attempt to determine the cau.se 
of the anomalies, we approximated the 
spectral reflectivity curves of the anom- 
alies and adjacent materials in the fol- 
lowing manner: The reflectivity for 
typical Mojave Desert materials (soil, 
playa, basalt) represented in the ERIM 
images was obtained from the literature 
lllovis, 1966); these spectral reflectivity 
values were assigned to the corresponding 
units in the image channels and u.sed to 
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Figure 3. Skykb S-I92 images of a part of track 6 in. southern California and 
Nevada: A, channel 7 (0.78 to 0.88 /tm). Black Cone located at 36°49' N., 116“33.5' 
W., neither Black Cone nor cinder cone Just soutliwcst of it show higli reflectivity; 

B, channel 11 (1.55 to 1.75 jum), note higli-refiectivity anomaly on botJi cones; 

C, channel 7 (0.78 to O.SS jum), compare with D; D, channel 1 1 (1.55 to 1.75 flm), 
notehi^i-reflectivity anomalies at Antelope Mountain (a), north of Antelope 
Mountain (b), nortli of Excelsior Mountains (c, d), and Black Cone (e). 
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calibrate the step wedges that accom- 
panied the images. Next, in order to 
assign spectral reflectivity values to the 
unknown material, a logarithmic inter- 
polation was made between the known 
values on the step wedge. Curves were 
then constructed for several of the 
anomalies (Fig. 4A). The accuracy of 
these curves is probably within 1 0 to 
20 percent. 

CONCLUSIONS 

There is little similarity between the 
reflectivities derived front the bRIM 
images (Fig. 4A) and spectral curves 
for common intermediate and mafic 
rocks (Hunt and others, 1971 ; 1973, 
written commun.). However, inspection 
of the reflectivity of the oxides (Hunt 
and others, 19731 hematite, iiinonitc, 
and goethite shows an increase in 
reflectivity in the near infrared. Three 
particle si?.es are presented in the spec- 
tral data of Hunt and others (1971 ; 

Figs. 4B, C, D), but it is only in the 
smallest size range (0 to 5 ^im) that the 
curves approach the very high reflectivity 
of the anomalies. 


That the mere presence of iron oxides, 
in any particle size, is not sufficient 
to produce the anomalies is further 
demonstrated in the ERIM scanner 
images of the brick-red Aztec Sandstone 
in the Spring Mountains (lat 35“55' N., 
long 1 15°2S' W., Fig. 1 A), which showed 
no higSi-reflectivity anomalies. 

A reasonable geologic explanation 
for the high-reflectivity anomalies has 
been proposed by botli Ray Wilcox and 
Richard White of the U.S, Geological 
Survey (1974, oral commun.). They 
suggested that the iron oxides (most 
likely hematite) were produced by late 
eruptive oxidation of mafic minerals 
and by precipitation of sublimate during 
posteruptive fumarolic steaming. This 
would account for the higli reflectivity, 
especially in view of the extremely fine 
particle size (probably <5 juin according 
to White). 

It is also interesting to note that the 
crater floors of tlie cones (Fig. 2) 
typically do not show the anomaly. This 
is consistent with ob.servations made by 
Wilcox, who noted that the sublimates 
are commonly absent in Hie floors of 
craters wliere early cooling, and hence 



plugging, of the plumbing systems im- 
mediately below the crater floors 
occurred prior to the posteruptive 
steaming. 

To date, 107 examples of the 
anomaly, as described above, have been 
observed on aircraft and satellite images 
and correlated with known geology. Of 
these anomalies, 39 occur in basalt 
(18 basalt flow’s, 21 basaltic cinders), 

63 in andesite, 2 in alluvium, and one 
each in intermediate-to-felsic rocks, 
felsic rocks, and limestone. This last 
occurrence is mapped as limestone on 
the geologic reconnaissance map; how- 
ever, the Skylab S-190-A photographic 
images show the area to be highly dis- 
sected and the anomaly to be reddish 
brown. 

Although these anomalous areas 
are at present only a geologic curiosity, 
it is likely that the technique used for 
discovering them can also be used for 
rapid discrimination of areas of post- 
eruptive steaming and sublimation, in 
addition, investigations are presently 
under way to .see if the intensity of the 
anomaly can be correlated to the age of 
the major geologic units. 
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Figure 4. Comparison of (A) spectral reflectivities generated from ERIM 
scanner and spectral reflectivities of (B) hematite. (C) gocthitc, and (D) limonite 
from laboratory data of Hunt and others (1971). 
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GEOTHERMAL RECOHNAISSANCE FROM QUANTITATIVE 

ANALYSIS OF THERMAL INFRARED IMAGES 

Kenneth Watson 

U.S. Geological Survey 
Denver, Colorado 80225 

ABSTRACT 

Remote sensing techniques may be useful In recon- 
naissance exploration of geothermal systems particularly 
because of our present limited knowledge of these systems. 

Although thermal infrared surveys provide a direct measure 
of surface heating, the variations In surface temperature 
due to geologlc-topograrhlc factors tend to mask the 
effects of geothermal heating. A thermal model is used 
to demonstrate that measurements of the mean diurnal 
temperature together with reflectance measurements and 
slope information can be used to map the geothermal flux. 

Theoretical calculations show that the error estimate 
of the mean diurnal temperature falls off rapidly 

with Increasing sample measurem.ents per diurnal cycle 
(650 HFU at 1 sample/cycle and 75 KPU at 3 sample/cycle). 

The thermal model was also used to examine the 
optimum times in the diurnal cycle In which to obtain an 
estimate of A thermal and a reflectance image of 

the Raft River, Idaho, area, acquired at the optimum 
time for 1 sample/cycle shows the presence of warm thermal 
anomalies which do not coincide with changes in the reflec- 
tance image. Probe measurements in the vicinity of one 
of these areas confirm the presence of such an anomaly. 

1. INTRODUCTION 

Geothermal energy is one of the alternate sources of power currently being 
evaluated as a supplement to conventional fuels. Because of increasing fuel short- 
ages, it is critical that reconnaissance techniques be developed which will enable 
us to understand geothermal systems so that we can evaluate their energy potential 
and develop extraction technology. At present, little is known of this potential 
and, until recently, exploration techniques were similar to those used in oil ex- 
ploration during the last century, which relied on location of oil seepage at tne 
surface (Banwell, 1970). 

Remote sensing techniques appear to have some promise for reconnaissance use 
in geothermal exploration. The work conducted in this study is directed at the use 
of thermal infrared surveys to detect geothermal heat and, in particular, the use 
of calibrated thermal images, acquired at optimum times during the diurnal cycle, 
to determine the geothermal detection threshold from thermal infrared surveys. 

2. CLASSIFICATION OF GEOTHERMAL AREAS 

The detection of geothermal anomalies by remote sensing will depend on the 
recognition of effects due to abnormal heat fluxes at or near the surface. White 
(1965) classified thermal "systems" into four broad categories: 1) normal systems 

[1.2 to 1.5 heat flow units (HFU)*; e.g., stable shield areas], 2) higher-than- 
normal systems [1.5 to 3-0 HFU; e.g., floor of the Gulf of California and its north- 

» 1 HFU H 1x10-0 cal/cm-Vsec 

1919 


ORIGINAL’ PAGE IS 
OP POOR QUALiry 





ward extension], 3) composite hydrothermal systems [3-0 to 15 HPU; e,g. , floor of 
the Salton Sea, California, and Larderello, Italy], and 4) hot spring areas [15 to 
1500 HPU; e.g.. Upper Geyser Basin, Yellowstone National Park, Wyoming; Walrakei, 
New Zealand; Steamboat Springs, Nevada; The Geysers, California; areas in 
Ahuachapan, El Salvador; and Reykjavik, Iceland]. 

Recognising the importance of separating the solar flux from the geothermal 
flux, Hase (1971) proposed three classes of geothermal anomalies with intervals 
based on the absorbed solar flux of 2000 HPU (Sellers, 19^7), snd his computed 
value of 240 HPU for the diurnal flux at a depth of one wavelength of the diurnal 
heating. 

A physically more realistic classification of anomalies for thermal surveys 
would recognize the differences among fluxes- detectable by Instrument, fluxes 
detectable on a single thermal image, and fluxes determinable from analyses of 
several images during a diurnal cycle. A detector thermal resolution of 0.5"C can 
be equated to a surface flux equivalent to about 50 HPU. The minimum heat flovJ 
anomalies (reported in the next section) determined from examination of individual 
thermal infrared images (generally obtained at night to minimize solar heating 
effects) have been 100 to 900 HFU’s (the lower value is optimistic, as will be 
shown later in otir analysis). 


3. PREVIOUS WORK 

Remote sensing techniques have been used in geothermal areas both for quasi— 
quantitative measurements of anomalous heat or flux and as a tool for supplementing 
conventional geologic and geophysical mapping. Those remote sensing techniques 
which appear applicable to the measurement of heat flov; include: 1) direct measure- 

ments of heat thermal infrared imagery and microv:ave radiometry, and 2) indirect 
heat measurements, i.-i., photography for snowfall calorimetry and mapping thermally 
stressed vegetation. 

Infrared radiometers and scanners have been used to map a wide variety of 
geothermal features over the last 15 years. A review of the literature for the 
years 1958-68 by Friedman and Williams (1968) indicates that a preponderance of this 
work concentrated on volcanic areas . The present paper uses a more limited defini- 
tion of "geothermal” vihich excludes active volcanism, because the detection and 
monitoring of these processes are treated more as hazards than as resources. 

One of the earliest direct measurements of surface temperatures in a geothermal 
area, using thermal IR (infrared), was done in 1961 at Yellowstone National Park by 
USACRREL (U.S. Army Cold Regions Research and Engineering Laboratory) and the 
University of Michigan. McLennan and Morgan (1965) reported that geothermal 
features could be clearly seen, but not identified by type, on the thermal infrared 
image obtained in this study. Yellowstone National Park v;as subsequently the site 
of a variety of geothermal investigations: Hiller (1965) used thermal IR in the 

identification of v/arm areas and recommended repetitive overflights in order to 
study dynamic changes in geothermal systems; VIhite and Hiller (I969) calculated 
that 150-500 HPU was the lo’wer limit of the detectable geothermal flux; Hase (1971) 
determined the detectability limit as 300-900 HPU. 

Strangvmy and Holmer (1966) detected a 6OO-HFO anomaly in the Lordsburg, 

N, Mex., area, using instruments capable, according to their estimation, of detect- 
ing a temperature difference equivalent to a 15— HPU anomaly. They also noted the 
importance of certain microclimatologlcal effects which produce anomalies on the 
radiometer scans; they suggested that microvjave observations, being less sensitive 
to such effects, v/ould provide better data. 

In Walrakei area of New Zealand, Hochstein and Dickinson (1970) estimated a 
detectability limit of 350 HPU for daytime imagery and 150 HPU for nighttime data. 
The factors vjhich they found affected the images Included solar radiation, wind, 
fog and condensation, and emisslvity. 

Palraason and others (1970) determined that the thermal infrared scanner was 
only moderately successful for the detection of anomalous heat flows in geothermal 
areas in Iceland. Their estimate of 300-700 HPU for the detectability limit closely 
approximates those determined for Yellowstone National Park and VJairakei. These 
authors stated that "with a calibrated Infrared scanner utilizing a DC restoration 
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circuit, and a known energy transfer function, by means [sic] of a finer amplitude 
slicing of the lew level anomalies, it may be possible to distinguish better between 
heretofore undetected low heat flow anomalies and the 'noise* caused by sources 
other than variations in heat flow to the surface." 

Gomez Valle and others (1970) reported that spatial resolution, rather than 
detectability, limited the thermal resolutions of the infrared images obtained in 
Michoacan, Mexico. Spatial resolution of f’ m caused the hot pools (v/ith a 1-m 
average diameter) to merge into their surroundings on the Images, therefore lowering 
the apparent temperatures. Improved scanner resolution and lov;er aircraft altitudes 
were suggested as soluclons. They also recommended a modification of tape readout 
levels to drop out all areas belov/ 30°C; such areas vrould be black on the film 
image. This modification vrould allov/ better discrimination of small changes in 
teinpex’ature above 30°C by "stretching" the temperature levels between SO^C and 
i00°C over the entire range of film sensitivity. The Michoacan report concluded 
x;hat thermal infrared studies should be combined vfith color infrared, electrical 
resistivity, heat flow, and hydrological and geoclieralcal studies in order to gain a 
better understanding of geothermal field. 

Moxham (19("9) reported that the abnormal radiance at The Geysers, California, 
is restricted to hydrothermal alteration zones and other local features and he 
found no evidence for a regional anomaly. He concluded that "The aerial IR 
technique of mapping terrain radiance may find its principal application in prelim- 
inary regional assessment of undeveloped geothermal areas." This suggestion proved 
Itself valid in the first small-scale geothermal reconnaissance survey undertaken 
in unexplored areas of Ethiopia and Kenya b,’ the United Nations Economic and Social 
Council, Committee on Natural Resources (1971). Flying a scanner vfith a minimum 
resolution of 6 in, at altitudes of approximately 3,000 m, the U.N. subcontractor 
found 105 anomalies of at least 10 °C in Ethiopia. Helicopter-borne field investi- 
gators found many of the anomalies to be hot springs, fuma xles, and warm fault 
traces. Similar features were detected in Kenya from a sco^iner with a 2-ra resolu- 
tion floim by the subcontractor at 915 m. 

Indirect methods for determination of anomalous heat flovi in geothermal areas 
Include color IR photography, snov;fall calorimetry, and vegetation mapping. ^ Color 
IR photographs were used by Dawson and Dickinson (1970) in Nevx Zealand for mapping 
heat losses in predominantly conductive parts of geothermal areas. They concluded 
that their method of contouring heat effects from color differences of vegetation 
on enlarged color IR images could shov; 2“C temperature differences at 15-cm depth. 
They estimated that the accuracy of the method vxas vrithin lOS of that of a ground 
survey and took only about 6 % of the time. 

Vfhite ( 1969 ) used snowfall calorimetry as a rapid means of determining both 
conductive and convective heat flow viithin a factor of tvio. V/hite mapped the Old 
Faithful and Sulphur Cauldron areas of Yellovfstone National Park using heat flovi 
contours ranging from f(50 to 5500 HFU. His method involved mapping of snow-melt 
lines combined vilth direct measuring of heat required to melt a known volume of 
snovx along the melt-line contour. Snow must accumulate to a depth of 10 to 15 cm 
during a 2 Ji-hour period for this technique to be practicable. 

Vegetation changes near the edges of the geothermal anomaly at Lordsburg, N. 
Mex., were identified on black and vihlte aerial photographs by Kintzinger (1956). 
These vegetation changes coT-responded to heat flow contours mapped from temperature 
probe data obtained an 1— m depth. Miller (1968) mapped similar stunted and dying 
vegetation around a geothermal area in Yellowstone National Park, but concluded 
that most of the vegetation changes viere caused by acid produced in the soil rather 
than directly by high heat flow. Perhaps both effects are important. 

4. THE THERMAL MODEL 

Various thermal models have been developed to analyze the surface and near- 
surface temperature behavior of the ground through a diurnal cycle. In a model 
developed by Jaeger (1953) for analysis of lunar surface heating a general Laplace 
transform solution v;as used which allowed him to treat nonlinear radiative transfer 
from the surface of a half space. This technique was applied (Hatson, 1971; VJatson 
and others, 1971) for Interpretation of nirberne infr„T’ed images, and later extended 
(Watson, 1973 ) to the case of periodic heating of a layer over the half space. The 
Laplace transform solution, however, is time consuming and requires iterative 



solutions; in addition, the physical significance of the terms is not clear. 
Therefore a linearised approximation to the radiative transfer law was introduced, 
the Fourier series approach was used (Lettau, 195^; Jaeger and Johnson, 1953), and 
key oases, thus analysed, were checked using the more precise Laplace transform 
solution. The numerical agreement was adequate for our purposes. 

The model assumptions are those previously introduced (Watson and others, 
1971): periodic, one-dimensional, and heat flow with constant thermal properties. 

The thermal behavior can thus be described using the diffusion equation. 

K at 


surface boundary conditions, <at x=o 
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v(x,t) is the ground temperature at distance x below the 

air-ground interface, at time t measured from local 
solar noon. 

K,k are, respectively, the thermal conductivity and thermal 
diffusivity of the ground 

e,A are, respectively, the average emissivlty and albedo of 
the ground surface 

a is the Stefan-Boltzraann constant 

Tgj^Y effective radiant temperature of the sky 

3 is the solar constant 
0 

f is a cloud cover factor to account for daily variations 
in the Insolation 

H(t) is the solar incident flux factor at time, t 
Q is the geothermal heat flux 


A linearized approximation to the radiative transfer term in equation (2) 
results in 
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Now the solar incident flux term is 
H = M{Z) cosZ’ 

where M(Z), the atmospheric transmission function, is approximated by 
M(Z) = 1 - 0.2 VsecZ 

Z is the zenith angle 
cosZ = cosi cos6 oosv;t + sini sind 
vihere 6 la the solar declinacion 
X is the latitude of the site 

w is the diurnal angular frequency 


and cosZ' = cosJ.' cos6 cos(wt)' + sinX’ sln6 


where Z' 1b 
X' = 
Cwt ) ' = 

d - 


the zenith angle for a nonflat surface 
X - d sin (<i> - ^) = X+d costf 
wt + d cos = wt + d siniji 

topographic slope 


- direction of slope 

The general method of solution to Eq. 1, subject to the boundary condition 
Eq, 2 and linearized radiative transfer (Eq. 3') is found in Carslaw and Jaeger 
(1959, p. 70-71). 


We are interested only in the surface temperature : 
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where the harmonic coefficients Cn and Dn are functions of the solar incident flux, 
H, the sky radiation^ Tojjy, the eraissivity, e, and the thermal properry ratio Kic"^ 
which is commonly called the thermal inertia. The mean diurnal temperature (Vijc) 
is obtained by integrating the ground temperature over the diurnal cycle, thus 
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Rearranging the terras yields the geothermal flux 
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The last terra on the right-hand side of equation (5) is a function of the 
ground albedo and slope and can be directly determined from quantitative reflectance 
measurements and topographic map data. If the surface bidirectional reflectance 
can be modeled satisfactorily according to the Lambert law then airborne reflectance 
measurements acquired at three separate times during the day can be used to compute 
this term. The three measurements are needed to solve for the unknowns: albedo, 

slope, and direction of slope. 

The TskY fcei*m can be directly determined using a calibrated sky radiometer. 

Its variations over the site should be negligible under clear sky conditions and a 
single value should suffice. The remaining unknown term in the equation is the 
surface eraissivity. It may be possible to make thermal measurements in tv?o spectral 
bands of the 8-l4 ^m windov: (Vincent and Thomson, 1972; Vincent and others, 1972) 
and use a spectral ratio to derive an empirical correction for eraissivity variations. 
In soil-covered areas this spec'tral contrast will be low and spectral correction 
may be Impossible. Fortunately, in these cases the departure from black-body 
emission Is a second-order effect and an eraissivity correction is unnecessary. 


5. ESTIMATION OF V,^^ (MEAN DIURNAL TEMPERATURE) 

The mean diurnal temperature (V,^c) be obtained by measuring the average 
ground temperature with a calibrated thermal scanner from repetitive aircraft 
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flights throughout the day and night. Theoretical calculations were made of the 
number of these measurements which are necessary to estimate V^c to various 
accuracies; the results are shown In figure 1, The temperature difference (AV) 
between this computed average temperature and the true was converted Into 
equivalent heat flovi -units (AQ) from equation (5): AQ = SeaTgj^y AV, 

The error bars mark the standard deviation in AQ for a 10° slope oriented 
in the four cardinal directions. The calculations were then repeated but with 
the sampling times offset from local solar noon by various time lags. The results 
are plotted in figure 2. 

It appears j from the data of figures 2 and 3j that optimum results should be 
obtainable from three or four samples per diurnal cycle; also, if data are acquired 
only once in the diurnal cycle then the optimum time is either sunset or shortly 
after sunrise. Figure 2 illustrates the large heat flux eri’or inherent in data 
acquired for geothermal analysis in the predawn. Actually, tvro types of errors 
are present: 1) an absolute error because the predawn temperature is not a good 

estimate of the average diurnal temperature, and 2) a relative error because the 
predawn temperature is still moderately dependent on thermal inertia (the de 
temperature is not). 


6, INITIAL APPLICATION OP MODEL 

Thermal emission and reflectance images of the Raft River, Idaho area v?ere 
acquired in early November of last year. Much of the area is flat valley alluvium 
(figure 3), covered with sparse scrub vegetation (figure which on oblique viev/ 
appears to completely cover the ground (figure 5)- Hot wells exist in the area 
and at one location (figure 6) the hot v/eli water flows out on the ground. There 
is, however, no altered, steaming ground in the area. 

The images (figure 7) were acquired at three different times in the diurnal 
cycle (08:39, 09:*J8, and 15:5“* MST) using the U.S. Geological Survey dual channel 
RS-14 a scanner in its first operational test. Optical alignment problems, since 
corrected, prevented a complete evaluation of the model; however, the data are 
satisfactory for illustration of the major points in the analysis. 

Data acquired at 09:lf8 MST most closely approximate the dc temperature. 
Comparison of the thermal channel (8-l4 um) data with the reflectance channel 
(0.5-1. 3 um) data (figure 8) shows the presence of warm thermal anomalies which 
do not coincide with the reflectance changes. An anomaly is particularly notice- 
able in the vicinity of the borrow pit where the warm area extends beyond the 
high reflectivity anomaly. Measurements with 1-m probes confirm the anomalous 
warmth of this area (Paul Williams, oral coramun., 1973). The image also defines 
another warm area near a gravel lens through which v/arm water rises to the surface. 
These hypotheses vxill be examined more intensively using heat flow measurements 
and repetitive thermal scanner flights. 
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Figure 1. - Equivalent heat flux error, AH versus the number of samples per 
diurnal cycle. The mean values (o) and the standard deviation (bars) were 
computed for a series of time lags from local noon. Values were computed 
for the circumstances of the Raft River: latitude 42*, solar declination 

-15”. 
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Figure 3. - Raft River area; view looking west. 



Figure iJ. “ Vertical view of grouml showing scrub 
vegetation and valley alluvium. 
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Figure 5. - Oblique view showing vegetation. 



Figure 6 . - View of hot water overflowing from well. 
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7. - Thermal Images of the Raft 
09:^8, and e) 15-‘5^ F.ST. Loca 



River area at a) 08:39, 
sunrise occurred at 07:00. 
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Pleure 8. - Thermal (a) and reflectance (b) images of the Raft River at 09 • MST, 
corresponding to one of the optimum times for geothermal mapping (see figure 
2) for one sample per cycle. The borrow pit (B) is an area of both high 
reflectivity and high temperature, indicating a region of anomalous reflec- 
tance and thermal inertia. The westward extension of the thermal anomaly 
beyond the anomalous reflectance region is consistent with increased geother- 
mal heating. 
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Geologic Applications of Thermal Infrared Images 
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/tto/ncf-Thenna] inrnired images provide inTormation about the 
near-surrace physical sute of geologic mateiials, particularly, the den- 
sity, water content, and heat transfer. Nontenestrial planetary studies, 
conducted at fairly coarse resolution, have been useful primarily in 
determining the distribution of rock fngmenis. Terrestrial studi'-s, 
conducted from satellite and aucraft at coarse to fine resolutions, have 
been successful for monitoring effusive volcanism, delineating areas of 
steaming altered ground and hot-spring activity, detecting fractures 
expressed hydrologically and topographically, and distinguisliing a 
variety of geologic materials with physical and compositional 
differences. 

Interpretation of thermal images is complicated by the various types 
of physical processes involved and commonly requires an assessment of 
many different factors. A simple theoretical model was used in this 
analysis to provide quantitative assessment of some of these facton, to 
predict optimum times to acquire thermal data, and to determine 
quantitative values of various properties of terrain. 

Two geologic applications were studied in some detail: geothermal 
mapping and thermal inertia mapping. Initial results indicate tliat both 
techniques have considerable potential, especially in reconnaissance 
studies. These data were acquired under optimal meteorological con- 
ditions and at sites where the geologic materials were well exposed. A 
realistic assessment of the limitations of these techniques must await 
future studies. 

Introduction 

O PTICAL-MECHANICAL scanners provide the means to 
monitor the temporal and spatial variations in the natu- 
ral thermal emission from planetary surfaces. These 
variations, observed at a great variety of altitudes from low- 
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flying aircraft to orbiting satellites, can be used to deduce 
physical properties of terrain materials and thus to infer 
geologic processes which have occurred. The interpretation 
of these data can range from direct visual examination of 
photographic recordings of the measured signal (using tech- 
niques developed in photogeology) to sophisticated digital 
processing using modeling analysis and pattern recognition 
techniques. These investigations are limited, liowever, by the 
complexity of the problem, in terms of both the physical 
phenomena and the number of different factors that influence 
the result. Of necessity, the interpretation of thermal infrared 
data has been based on quite simple theoretical models, in- 
volving very limiting assumptions and fairly ideal circum- 
stances (both meteorologically and geologically). Future re- 
finements in observational technique and theoretical analysis 
will be required in order to apply the technique routinely to 
the diverse terrain conditions found in nature. 

This paper describes the development of one theoretical 
model for analysis of the surface temperature distribution. 
The results provide a quantitative estimate of the effects of 
various geologic, meteorologic, and topographic factors and 
hence can be used as a basis for direct interpretation of the 
thermal images. In addition, the model curves can be used 
to predict optimum times to acquire the data in order to 
enhance specific effects and features. Finally, the model can 
be used to determine quantitative values of various properties 
of the terrain. This last aspect of the model has greatest 
impact in future geologic applications and receives the most 
attention in this paper. 
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The Ti£ERMal Model 

A mathematical model of the near-surface conductive heat 
transfer was used to develop a means to interpret thermal 
infrared images in terms of the relevant physical properties 
and processes. For simplicity, I assumed one-dimensionat 
periodic heating of a uniform half-space (a region bounded 
by a plane on its upper side and extending downward to 
infinity) of constant thermal properties. The temperature 
obeys the diffusion equation: 


du 

^ dx^ dt 


CD 


where v = u(x, f) is the temperature at distance x below the 
surface and time f, from local solar noon, and k is the thermal 
diffusivity of the half-space. The solution to C I ) for periodic 
heating of angular frequency w is [4, p. 8 1 ] 


uCx, f) = Z) exp (~k^/nx) cos {nut - “ k-^x) 12) 

n=0 

where and are arbitrary coefficients and fc=Vt3/2K is 
the wavenumber of the first harmonic. The arbitrary coeffi- 
cients are evaluated by expressing the surface boundary con- 
dition as an energy balance between incoming solar and sky 
radiation, outgoing ground radiation, and conduction into 
the ground. For the moment, I shall ignore other heat trans- 
fer mechanisms (atmospheric conduction and convection, 
convective heat transfer in the ground, and latent heat effects 
due to evaporation or condensation of water). 

Thus at X - 0 


-K 


9u(Q, t) 
dx 


= -eau* 


+ / 


(3) 


where I is the absorbed incoming radiation composed of both 
the short-wavelength solar flux (less than 4 jam) and the long- 
wavelength atmospheric flux (greater than 4 jum), e is the 
mean emissivity of the surface, a is the Stefan-Boltzmann 
constant, and K is the thermal conductivity of the half-space. 
The 'quadratic temperature term precludes direct substitu- 
tion of (2) in (3). Jaeger [15] used Laplace transforms to 
develop a relation between the surface flu.x and tlie surface 
temperature, and solved for surface temperature by iteration. 


p m 

fii " ~ 7 = Us^rt-j+i> n~!,2, (4) 

Vrrr^i 

where /ff is the average flux into the ground in the nth interval, 
P is the thermal inertia of the ground (P = iI/'»yK"), r is the 
period of the heating flux, is tlie average surface tempera- 
ture in the rth interval, and the <p's are a set of numerical 
coefficients determined solely by tn, the total number of in- 
tervals in r. The term /„ which is equal to the right-hand 
side of (3) contains a term thus requiring an iterative solu- 
tion for tij. Tliis approach has been used both to analyze 
planetary data [161 and to interpret aircraft thermal images 
[44] , [45] . The technique has been e.xtended to tlie case of a 
layer over the half-space [43] in order to examine the effects 
of surface coatings and lichen cover. 

A second approach to Uie nonlinear transfer problem is the 
method of finite differences [46]. Care must be taken in 
selecting the spatial and time increments to ensure that con- 
vergence is achieved. Both the finite difference and the 


Laplace transform technique obscure the physical significance 
of the terms and can be computationally time-consuming. It 
was decided to linearize the radiative flux term in the bound- 
ary condition and then check the numerical results with the 
more exact Laplace transform solutioir. The comparison was 
satisfactory for terrestrial diurnal temperature ranges. The 
linearized solution [9] is too inaccurate for the temperature 
variations encountered on the lunar surface [ 16] . 

The incoming radiation term / is composed of tjie short- 
wavelength radiation from the sun {modulated by the at- 
mospheric transmission) and the long-wavelength radiation 
from tile sky. This latter term can be approximated by an 
expression urs^y where Tjfcy is an effective sky radiance 
temperature; hence the absorbed flux is eoT^y. The long- 
wavelength terms of the right-hand side of (3) can then be 
linearized as 

eov^ - ~ deaPs^y (w Psky), for — « !, (5) 

^sky 


The solution of the diffusion equation (1) subject to the 
boundary condition (3) modified by (5) can be obtained by a 
simple substitution [4, p. 70] ; let 


where 


Then 


0{x, f) = 0“ 
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subject to the boundary condition at x = 0 

0 = 7’sfcy+7:r/A'h (8) 

where 4 is the short-wavelength absorbed flux. This term is 
a function of the ground albedo A, the soiar declination 5, 
the site latitude A, and the local slope and can be e.xpressed 
as follows: 

4 = (1 - ^) SoCMCZ) cos t (9) 

where So is the solar constant, C is a factor to account for the 
reduction in solar flux due to cloud cover, M{Z) is the atmo- 
spheric transmission as a function of tlie zenith angle Z, and 
Z* is the local zenith angle for an inclined surface. 

The atmospheric attenuation can be approximated by a 
square root of zenith angle law [ 1 ] ; let 

itf(Z)= 1- 0.2Vi^ (10) 

where 

cos Z = cos X cos 6 cos tor + sin A sin 5. (11) 

The local zenith angle z’ can be computed from 
cos Z' = cos d cos Z - sin d { sin 0 cos 5 sin tot 

- cos 0 sin 9 cos A - sin 5 sin A cos cjt} (12) 


where d is the surface slope angle measured downward from 
the horizontal, and 0 is the azimuth of the slope angle mea- 
sured counterclockwise from north, For convenience, I shall 
define an additional parameter H(t) which is the local 
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Fig. I. Diurnal temperature curves for varying (a) thermal Inertia (in cal/ctnVs'''^), (b) albedo (in fractions <1), (c) geothermal heat flux 
(in HFU), and (d) emissivity (in fractions <I). Times are local solar values. The mean diurnal temperature DC is indicated by a hori- 
zontal line for each curve. Note that this value is constant for varying thermal inertia [see 1(a)] and lliat the amplitude of each curve 
is reasonably independent of geothermal flux [ see 1 (c) ] and emissivity [see 1(d)]. The vertical arrows near the time axis mark the sun- 
set and sunrise times. Fixed parameter values used in the computation of these curves are: inertia 0,03 cal/cm^s*'^. albedo 0.3, emis- 
sivity 1.0, latitude 30°, solar declination 0°, sky radiant temperature 260°, dip 0°, cloud cover 0.2, 


insolation; 

tMiZ) cos Z\ 

where tg and are the local sunset and sunrise times (mea- 
sured from local solar noon), and are chosen so that for 
“fjj > f > fg, cosZ> 0, and cosZ*> 0. 

Hence 

= (13) 

and the boundary condition at JC = 0 (S) becomes 

= Tsky + U - A)SoCfmh. (14) 

The solution to (7) subject to (14) is 


0C-C, 0 = Tsky 


(I-/l)SoC 

Kh 


♦ ^ exp i-k\/nx) cos (liCJf - - k\/nx) (15) 

where and e„ are the amplitude and phase of the harmonic 
components of the local insolation H. 


The surface temperature can then be determined from equa- 
tion (6) by integration [4, p, 70] 


u(0, f) = h 



0(i7, r) exp (-/JT?) dn 


~ ^sky 


(1- vDSpC 

K 


z 


A„ cos (nojt- e„ " Sj,) 
^/{b + ky/n)~ + (k.yjnf 


(16) 


where S„ = tan“‘ {k\/nlh + [ky/n )). 

The effect of a geothermal heat flux Q can now be intro- 
duced by adding a second solution u = Qx/K + Q/Kh which 
satisfies the boundary condition (3) and the differential equa- 
tion (1). 

Recall tit it 

k =\/oii2K 
h = 4ecrrsky/A:, 


Introduce r = P\Ji^, where t = 27r/to and s = UK, then 
kyjnlh = ryjnls and 

yjkn + 7l^ = hyJ(jry/n(s)‘^ + 1. 
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Hence 


u(0, f) = Ts^y + Q/s + 0 - A)SoC 


where 


■ 

r» =fl 


cos (ncjf ~ ~ Sft) 

y/(s+ rv^)^ + iry/n)^ 


(17) 


5„ = tan"^ (fV^J + i-y/n) 
r = Py/irJr 
s = 4effrsly. 


Fig. 1 illustrates the diurnal temperature behavior exhibited 
for variations in four important geologic parameters; thermal 
inertia, albedo, geothermal flux, and emissivity. 

The mean diurnal temperature Udc “ computed by inte- 
grating (17) over the diurnal cycle 



Q 

u(0, 0 dt = 7*sky + — + (1 - A)SoCAq cos 6q/s 


(18) 


where Ao cos Gq = (1/r) /J" Hit) dt. 

it is important to note that ujc >s independent of the ther- 
mal inertia of the ground and can be used in conjunction 
with albedo measurements and topographic information to 
derive the geothermal flux Q. 

The day-night temperature difference (Ay) is computed by 
forming the difference between the noon and midnight 
temperatures. 


Au = u(0, 0) - y(0, r/2) 


= 2(l-^)SoC ^ 


A^m-l C0s(G2„,-i +52;„-i) 


= i s/is -^rs/lm - 1)^ + {ry/2m -T)^ 

(19) 


The quantity Au/(i ~ A) which is a function of the thermal 
inertia P (through dependence on r) can be determined ob- 
servationally and thus used to map thermal inertia variations. 


GEOTHERMAL MAPPING 

Thermal images display the spatial surface temperature 
variations of the ground and thus thermal infrared mapping 
is an obvious technique to consider for reconnaissance explo- 
ration of geothermal areas. Success has been achieved in 
monitoring worldwide areas of effusive volcanism (for exam- 
ple, in Hawaii [10], Iceland [25], and Italy [5]) and in 
delineating areas of steaming altered ground and hot-spring 
activity in the western United States (for example, Yellow- 
stone National Park (191, [20], The Geysers, Calif. [21]). 
Considerable reservations remain, however, about the general 
application of thermal mapping to the detection of lower level 
geothermal anomalies principally because of the noise effect 
of gcologic-hydrologic-topographic variations [2j. Estimates 
which have been made of tiie minimum detectable geothermal 
flux are based largely on the assumption that predawn mea- 
surements are optimum, because they minimize the effects 
of solar heating: Hochstein and Dickinson [13]. 150-350 
HFU;‘ Wliile and Miller [47], 150-500 HFU; Hase [12], 


* I HFU= IX 10“« cal/cm*/s. 


^ CasD 0 



THERMAL ALBEDO EMISSIVITY SLOPE 

INERTIA 


Fig. a. Equivalent heat flux error SQ computed from temperature dif- 
ference values between standard and comparative values for different 
geologic-topographic properties: thermal inertia P, albedo A, eniis- 
sivity e, and slope d. Tlie assumed standard values are P~ 0.03 
031 /^ 75 "*, A - 0.3, e= 1.0, and d = 0° and comparative values are 
P = 0,04, A = 0.2, e = 0.9S, and d = 10°. Case a is computed from 
differences in the predawn values arid ,‘s a measure of the “noise” 
present on predawn images which are used to make relative estimates 
oi the geothermal flux. Case b is computed from differences between 
predawn comparative and mean diurnal DC standard values and is a 
measure of the absolute error in predawn measurements. Case c is 
computed as in case b but using midmorning (10 A.M.) comparative 
values. Case c illustrates the signincant reduction in the "noise” ef- 
fects that can be achieved by correct selection of data acquisition 
times. 


>240 HFU (theoretical) and 900 HFU (measured); and 
Palmason et al [25], 200-700 HFU. It is instructive to 
compare these values with quantitative estimates of the ef- 
fects of various geologic-topographic factors. The compari- 
son is made more meaningful by converting temperature dif- 
ferences into equivalent heat fluxes, and this is done by a 
theoretical computation of the temperature change which 
occurs for variations in the geothermal flux. The results 
(Fig. 2) illustrate that natural variations in geologic properties 
or topography can easily overwhelm internal heat fluxes of 
several hundred HFU in all but the flattest and most homo- 
geneous terrains. 

The next step is to see whether some analysis scheme can 
be developed which will reduce the effect of some of these 
factors. In the section on the thermal model (18), 1 demon- 
strated a relationship between the mean diurnal temperature 
ude internal heat flux (0), Ude is independent of 

the thermal inertia of the ground and the equation provides 
a formulation to correct quantitatively for variations in albedo 
and slope. To compute Q there remains the problem of de- 
termining yjc. albedo, and emissivity by remote measurement. 
A computation of Vde can be made by measuring the ground 
temperature at several times in the diurnal cycle. Tlie number 
of times required can be evaiuated by comparing the mean 
of several measurements with Ude and converting this tem- 
perature difference into the equivalent heat flux. 
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Let 




1 

Si;(m) = “ V o( 0 , ffc)- Ode- 


The calculation of fiu can be further refined by the introduc- 
tion of a variable lag ( 0 ) measured from local noon. 

Let f* = 0 + (fc - 1) rfm, where 9 < r/m. Therefore, 


oo 


5u{m,0) = a“^)5oCX; ^pm 

p=0 


cos (pm CO 0 - Ep„, ~ dpf„) 
-\/{s +r\/pm)^ + (r-\/pm)^ 


( 20 ) 


Now from (17) it is possible to convert 6 u into an equivalent 
heat flux (SQ), he*, 5Q ==s5u. The results, plotted in Fig. 3, 
show that the estimate of ude can be significantly improved 
by increasing the number of samples to 3 or 4 samples per 
cygle. The gain achieved by further increase to 5 or 6 samples 
per cycle is largely illusory, because measurement errors and 
calibration problems will most certainly exceed the gain. 

The ground albedo can be determined by combining reflec- 
tance measurements v/ith topographic map information, pro- 
vided that we have an adequate model of the bidirectional 
reflectance. For simplicity, we can assume that bare and 
sparsely vegetated soils can be approximated as Lambert re- 
flectors [ 8 ]. In that case, reflectance measurements at three 
times in the diurnal cycle can be used to correct for albedo 
and slope. 

Let the measured reflectivity at time tj be rj. Then 
r/ = v4cosZ/, / = 1 , 2 , 3 
where cos Zj is defined as in (12). Then 
A =r 2 hosZi 

where 

d =-tan“* {(xj cos A cos 6 + Zi sin 6 sin 6 )/ 

* (xi sin A sin S “ y 1 sin ^ cos 5 + Zi cos 6 sin 5 cos A)} 

^ = 7 T + tan"* {DF- E-sjD'^ F^(EF + DyjD'^ +E^ - F^} 
and 

£» = cos 6 {(.vay 1 ~ j:iy 3 ) cos A cos 5 

+ (yiZs “ yaZi) sin A sin 5} 

£■ = sin S cos 5 cos" A t^iZa - X 3 Z 1 } 

F = "Sin^ 6 sin^ 5 (XiZa - x^Zi) 

X{ = Tf cos C0f2 - cos COff 

y,- = r'l sin o)t 2 ~ sin cof,- 
ri = rilr 2 , i=l,3. 

The remaining unknown in the geothermal equation (18) is 
the emissivity of the surface. Actually, emissivity variations 
enter the problem in two ways: emitted radiation and kinetic 
temperature. The thermal scanner measures the effective 
radiative temperature of ground and this can be approximated 
by e*^**u, where u is the ground kinetic temperature, and e is 
the surface emissivity. 

But tiie ground kinetic temperature and hence its mean 
Ude are also functions of emissivity. Because the kinetic 
temperature decreases with increasing emissivity while the 
radiative temperature increases, the two efi'ects work against 
each other and reduce the net effect of emissivity variations. 

^ a 



NO. OF SAMPLES 

Fig. 3, The equivalent heat flux error (in HFU) versus number of 
samptes/diurnal cycle. The mean and standard deviations are com- 
puted over variable time lags (computed from local noon). 



Fig. 4. The average emissivity (8-14 pm) versus the emissivity ratio 
(computed over 1-fim bandwidths centered at 10,5 pm and ll.s pm) 
for various igneous rocks whose spectral reflectance was measured in 
the laboratory (40 ]. The line is the computed quadratic least squares 
fit to the data; the standard deviation of this fit is 0.01 or approxi- 
mately I percent of the average emissivity. 


For most soils, the effect of fine grain size is to reduce the 
spectral contrast and raise the emissivity in the 8-I4-jim re- 
gion [36], producing a fairly constant value of 0.95. The 
most significant departure from this value occurs for silicate- 
rich soils or rocks, owing to the rcststrahlcn features near 
10 pm [7], [14], [18]. In these cases it may be possible 
to correct for emissivity variations by measuring the spectral 
emission ratio in two adjacent spectral bands. Spectral band 
ratios have been previously applied to enhance the spectral 
contrast of rocks and soils [37]-[391 , Using published labora- 
tory reflectance data of natural surfaces [40], 1 have com- 
puted the average 8-14-jLtm emissivity and the emissivity ratio 
in two adjacent bands. A quadrative least squares fit to these 
data (Fig, 4) shows that by this ratio the net emissivity can 
be estimated to an accuracy of I percent. Further examina- 
tion of natural surfaces on a scale comparable to aircraft 
-2. 
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Fig. 5. Tliermal (S-14 f«n [upperl) and reflectance (0.3-1.3 [lowerlj *mage pair for the Raft 
River, Ida., area. Ttie ini.ige acquisition time (09; 48) corresponds loan optimum time to estimate 
tlic DC temperature. Note the warm area in the vicinity of a borrow pit (just above center of upper 
image). It coincides, in part, with a liigli-reflectance anomaly tlower image] but extends beyond it. 


scanner resolution of several meters will be necessary to 
evaluate the practicality of this approach. 

An initial test of our tcchnictue to use for mapping 
geothermal heat flux vv.is performed at Raft River, Idaho 
[421. A thermal image, acquired at the optimum time for 
one sample per cycle, is shown in conjunction with the reflec- 
tance image obtained at the same lime (Fig. 5}. 

An interesting tlicrmal anomaly is present in the vicinity 
of a borrow pit. Part of the warm ground corre.sponds with 
a high-reflectance anomaiy associated with vegetation differ- 
ences. Ttic remaining anomaly lies outside this are.i and is 
not due to reflectance or slope effects. Because llie data 
we.rc acquired at an optimum time for geothermal mapping, 

^ 3-4 


the ground temperature is not strongly dependent on thermal 
inertia differences. It seems reasonable to attribute the en- 
tire anomaly to geothermal heating. This conclusion is con- 
sistent with probe measuremcnt.s at the site. 

It is tempting but still premature to speculate on the mini- 
mum detectable heat flux that can be achieved with this 
technique. The model results suggest that uje he esti- 
mated by repeated scanner overflights to an accuracy of .i 
few tens of IIFU. Conventional thermal scanners can coUect 
calibrated high-resolution images witJi an NETD (noise equiv- 
alent temperature difference) of 0.l“C which translates to 
less than 10 HFU, But this must be tempered by a comparison 
with other energy sources which have been neglected in our 
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analysis. For example, although we define the incident solar 
flux as short-wavelength radiation (occurring below 4 Jim), 
about 1 percent of the solar flux occurs beyond 4 jim and 
this amounts to about 300 IIFU. Under stable atmospheric 
conditions, radiation is the dominant heat transfer mechanism 
from the surface. At 300 K this flux is approximately 10 000 
HFU. Convective heat transfer, under moderately windy 
conditions, can equal this, and thus even under stable condi- 
tions, subtle microclimatologic and surface roughness varia- 
tions could produce heat flux errors of many tens of HFU. 
Then, too, errors in estimating the surface albedo, atmospheric 
transmission, and scattered sky radiation are significant be- 
cause these factors mor ulaie the solar flux. A realistic esti- 
mate will have to await more observational tests and a de- 
tailed comparison with near-surface probe data of the 
geothermal flux. 

Thermal Inertia Mapping 

Previous investigators recognized the utility of thermal 
infrared surveys to discriminate geologic materials based on 
differences in their thermal properties. The earliest geologic 
application was the telescopic observations of the moon [ 22 ], 
[ 26 ], [30], [33], (34), which demonstrated the presence of 
the lunar dust layer and the many anomalies caused by rock 
fragments in the vicinity of bright-ray craters. 

Terrestrial studies have been conducted both from aircraft 
and from satellite. Aircraft overflight data have been used 
to distinguish various geologic materials; olvine, basalt, and 
lapilli (ash beach deposits) and lacustrine carbonate deposits 
in the Mono Lake area, Calif. [11]; sedimentary bedrock 
versus windblown sand cover and individual siltstone and 
sandstone strata in the Imperial Valley, Calif. [31]; alluvial 
soils in the Rio Grande Valley, Tex. [23]; dolomites and 
limestones including facies changes between them in Mill 
Creek, Okla. [29] ; soils of varying density in the Caliente 
and Temblor Ranges, Calif. [48] . 

Terrestrial satellite observations have indicated thermal prop- 
erty differences which are measurable from space [24] , and 
which have been used to map different geologic materials 
[28]. Planetary probe data acquired of Mars [17] and Mer- 
cury [6] show the presence of some thermal property varia- 
tions at a fairly coarse resolution (50-100 km) which are prob- 
ably due, as in the Moon case, to surface and near-surface rock 
fragments. 

The surface temperature variation of a periodically heated 
homogeneous half-space is dependent on a single thermal prop- 
erty called the thermal inertia (4) and (17). This property, 
which is the ratio of thermal conductivity to the square root 
of the thermal diffusivity, correlates somewhat linearly with 
the density for most dry geolopjc materials (Fig. 6). In addi- 
tion, minerals and rocks high in silica (e.g., quartz, rhyolite, 
granite) tic above the correlation line (having somewhat higher 
thermal inertias than expected for their densities), and minerals 
and rocks low in silica and high in iron (e.g., olivini.. basalt, 
gabbro,dunite)lie below the correlation line (having somewhat 
lower thermal inertias than expected for their densities). Due 
primarily to their higher conductivity and specific heat, dolo- 
mites generally have iiigher thermal inertias than limestones. 
The addition of moisture to a dry soil results in a significantly 
larger increase in thermal inertia than should be expected for 
the increased density. This result is due to the replacement of 
air with water, significantly changing both the specific heat and 
the conductivity. For e.xample, a dry .sandy soil (density 
1.65 g/cm^, specific heat 0.19 cal/g, and conductivity 0.00063 



Fig. 6. Thermal inertia (cal/cm’/s*'*) versus density (g/cm’) for a 
variety of rock-forming minerals, rocks, and soils. The linear correla- 
tion line, bracketed by lines showing the standard deviation, shows 
the increase of thermal inertia with increasing density. Rocks and 
rock-forming minerals which are high in silica tend to fall above the 
line (see quartz). Those which are low in silica and high in iron oxide 
tend to fail below the line (see olivine). 

cal/cm®/cm/s/degree [4]) has a thermal inertia of 0.014 cal/ 
cm^/s*'- . On the other hand, a thermal inertia value of 0.043 
for a water-saturated soil (density 1.93 g/cm^) can be derived 
from the expression in [45] based on Sugawara and Yoshi- 
zawa’s equation for the thermal conductivity of porous mate- 
rial [35 ] . The saturated soil has a thermal inertia only slightly 
less than that of a sandstone-0.056 (density 2.3 g/cm^ ). 

In their natural state, geologic materials are commonly oxi- 
dized, stained, lichen covered, and covered with residual soils 
and various types of vegetation. Discrimination of geologic 
units by the thermal inertia technique depends on the extent 
to which these surface effects provide distinguishable thermal 
characteristics. A model consisting of an insulating layer over 
a homogeneous half-space was constructed [43] in order to 
examine the extent of some of these surface effects. Both a 
0.01 -cm "lichen” layer and a 0.1 -cm “soil” layer have a negli- 
gible effect; a 1-cm lichen layer and a 10-cm soil layer produce 
a surface temperature similar to an infinitely thick soil cover. 
Clearly then, the thermal inertia mapping technique is only ap- 
plicable to rock-type discrimination in areas of good geologic 
exposure such as those found in arid and semiarid regions. 

The technique has been applied using data acquired from the 
Nimbus meteorological satellites (8-km resolution) for part of 
Oman in the Arabian peninsula [ 28 ] . The test site was selected 
mainly because of exceUent exposures of diverse rock and soil 
types of large areal extent. The resulting map showed gross 
agreement with the major units on a reconnaissance geologic 
map. Some anomalies which did not corresp' nd to features on 
the reconnaissance map, subsequently, were found to match 
features on a more detailed geologic map or in a more detailed 
image obtained by the F.RTS spacecraft [27]. The technique 
has now been computer-automated and refined using atmo- 
spheric correction models. A preliminary map of the eastern 
part of the Arabian peninsula is shown in Fig. 7. Higher reso- 
lution data from existing and future satellites should make it 
possible to apply the thermal inertia technique to discrimina- 
tion of geologic materials and to study soil moisture conditions. 

Additional Afplications 

So far 1 have discussed those geologic applications which are 
based on a quantitative analysis of repetitively acquired images 
using a theoretical model. Other applications have been dcvel- 
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Fig. 7. Thermal inertia map of pari of the eastern Arabian peninsula from 20 to 27® N latitude and 
from SS to 60® E longitude. Contour values are multiplied by O.OOS to yield thermal inertia in 
cal/cm*/s‘'*. 


oped from visual examination of the image data. Reference 
has already been made to their use to monitor effusive voka- 
nism and to detect steaming ground and hot springs. Grouiid* 
water migration to the surface can cause, in addition U the 
obvious thermal property change, surface-temperature cooui'ie 
due to evaporation (from bare ground) or evapotranspiration 
(from vegetation). In some instances this enables recognition 
of fractures and fault zones [ 29 1 , ( 32 ] , (4 1 ] . In other cases, 
the surface cooling can lead to the recognition of areas of po- 
tential la.ad slippage (49), (3| . 

. - 



Thermal images showing contrasts due to differential warm- 
ing of oppositely facing surfaces commonly provide an en- 
hanced display of various topographic features. In some in- 
str.nces, this has resulted in the detection of lineaments which 
can provide information on the regional structure, and in some 
cases, may be related to the deposition of economically valu- 
able minerals (49). Topographic effects can also produce a 
thermal image pattern which provides geomorphic information 
ultimately related to such matters as rock type, fracturing, and 
erosion. In some cases, this prompted the conclusion that post- 
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sunrise images resemble low-sun angle photography [48]. In 
other cases, the images have aided in recognizing stratigraphic 
and structural detail [29] , (31) . Rowan et al. report enhance- 
ment of small-scale bedding detail in morning images and of 
alternating formations that form hogbacks and valleys in after- 
noon images. It may be possible to apply pattern recognition 
criteria to these images in order to develop discrimination cri- 
teria for different geologic materials, and to quantify the ef- 
fects of various geomorphic processes which act to modify the 
surfaces. 

Summary 

Thermal infrared images have been used to monitor effusive 
volcanism, delineate areas of steaming altered ground and hot- 
spring activity, detect fractures, and discriminate a variety of 
geologic materials. Because the interpretation of these images 
is complicated by the effects of many geologic, meteorologic, 
and hydrologic factors, a theoretical model was developed 
which can be used to predict optimum times for data acquisi- 
tion to determine quantitative values of various terrain 
properties. 

The model assumes one-dimensional periodic heating of a 
uniform half-space conductive heat transfer in the ground, and 
radiative transfer in the atmosphere. The resulting diurnal 
temperature variation is expressed as a function of geologic 
properties (albedo, thermal inertia, emissivity, and geothermal 
flux), topography (slope and azimuth),' and solar-atmospheric 
properties (transmission, sky radiance, site latitude, solar decli- 
nation, and cloud cover). Two important properties of the 
diurnal temperature variation are computed: the mean diurnal 
temperature and the day-night temperature difference. The 
mean diurnal temperature is independent of the thermal inertia 
of the ground, and can be used in conjunction with albedo 
measurements and topographic data to derive the geothermal 
flux. The ratio of day-night temperature difference to the 
complement of albedo is dependent, to first order, only on the 
thermal inertia and hence the ratio can be used to discriminate 
geologic material. 

Geothermal application of thermal images has been limited 
in the past by the geologic-topographic “noise” factors which 
can easily overwhelm all but the most intense anomalies. A 
theoretical analysis suggests that both thermal and reflectance 
image data should be acquired at least three times in the di- 
urnal cycle. A limited application of the model to data ac- 
quired at an optimum time illustrated the presence of a weak 
thermal anomaly wliich was confirmed by ground measure- 
ments. It is still premature, however, to estimate a minimum 
detectable flux using this technique. 

Thermal inertia mapping from Nimbus satellite data has 
already demonstrated the feasibility of discriminating a variety 
of geologic materials. Planetary probe data display thermal- 
property variations which are probably due to near-surface 
rock fragments. The dependence of thermal inertia on density, 
on moisture content, and, to some degree, on composition 
suggests that the technique should be useful to discriminate a 
great variety of geologic materials and conditions. Other geo- 
logic applications of thermal mapping include ground-water 
monitoring and geomorphic and stiuctural studies. 
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The microwave brishtness of heterogeneous materials is decreased by internal scattering. 
This study, based upon radiative transfer theory, shows that for low loss media darkening 
can be many tens of degrees, that the effect is greater for materials of lower dielectric con- 
stant, that darkening depends upon the dominance of scattering over absorption, that volume 
scattering tends to mask a temperature gradient, that darkening is slightly greater for V- than 
for ^f-polarized radiation, and that darkening increases slightly with view angle. For isotropic 
scattering, the analysis yields the brightness temperature as a function of the scattering albedo, 
the complex dielectric constant, thermal temperature, and the temperature gradient. 


INTRODUCTION 

Thermal radiation from rock, soil, or snow may be 
used remotely to infer physical properties such as 
temperature, dielectric constant, moisture content, 
or thermal inertia. Interpreting thermal-infrared 
emission (8 to 14 /t] is eased by assuming that the 
radiation is diHu.se and that it originates at the opti- 
cal surface. Interpreting microwave emission (10"^ 
to 10® m) is, in this sense, less simple. Surfaces are 
seldom sufficiently rough to produce diffuse radia- 
tion nor are they quasi-specular. Much of the radiated 
energy may originate several wavelengths below the 
optical surface forcing consideration of the physical 
characteristics of an emitting volume as well as the 
roughness of the surface. Pebbles in soil, ice glands 
in snow, vesicles in lava, stratification, or foliation 
scatters the radiation internally and so the emitter 
cannot be modeled as a homogeneous halfspace 
bounded by a rough interface. The contributing ef- 
fects of surface roughness, volume scattering, and 
layering must be unraveled before thermal tempera- 
ture or dielectric properties can be ferreted from 
thermal microwave emission. 

Siogryn [1970] soived the theoretical problem of 
emission from a uniformly layered halfspace, and 
Blinn et ah [1972] have observed the expected inter- 
ference patterns. Edgerton et al. [1971] discussed 
emission from a non-uniformly layered halfspacc in 
an attempt to explain the anomalous microwave 
signature of a metamorphosed, layered snowpack. 
Emission through a rough surface might be treated 
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similarly to problems'of radar backscatter [e.g., Beck- 
man and Spizzichino, 1963; Peake, 1959; Stogryn, 
1967]. 

This paper concerns emission from a halfspace 
containing randomly distributed scatterens embedded 
in a low loss dielectric and considers the relative in- 
fluences of a temperature gradient, the dielectric 
constant, and absorption. 

THEORY 

The model in Figure 1 is of an isotropic, low loss 
dielectric occupying the negative halfspace bounded 
at a = 0. The surface is quasi-specular (effectively 
specular at the wavelength of concern), the tem- 
perature within the dielectric is a function of depth 
(—a), and the complex dielectric constant is e* 
(= e' — u"). The dielectric contains randomly dis- 
tributed scattering centers of cross section o-, and 
density N. These are point scatterers and neither 
absorb nor emit but merely scatter. The radiation 
intensity at a point within the dielectric is specified 
completely by four parameters — those typically 
chosen are the Stokes parameters, Becau.se the radia- 
tion intensity is rotationaliy symmetrical about the 
z axis, phase information need not be preserved and 
the required number of Stokes parameters is reduced 
by half. That is, the total intensity can be written 
as tlte sum of horizontally and vertically polarized 
radiation, 4- 1%-. nie radiation intensity at wave- 
length Ao within the dielectric is an explicit function 
only of z, and of angle between the direction of 
energy flow and 4-z. If is the cosine of the angle 
then the intensity can be written Ipiz, fi) where 


3 . '^4 
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Fig. 1. Model for thermal emission for an isotropic 
dielectric occupying the halfspace z < 0. At a specific 
wavelength, /,(r. m) is the radiation intensity within the 
dielectric where r is the optical depth below the surface, 
and n is cos cosine of the angle between the direction 
of propagation and vertical. The radiation intensity above 
the surface is denoted /,(+ 0, where is view angle. 
The subscript p refers to either horizontally or vertically 
polarized radiation. 

polarization p is either H or V. The explicit depend- 
ence upon wavelength is suppressed. Note that <f>' 
and fi, whether p is primed or not, refer to direction 
within the dielectric. By convention, tj, is reserved 
for view angle above the dielectric. 

The equation for radiative transfer at z < 0 be- 
comes 

3/,(2, n)/ds = — a/,(r, p) + £,(r) + (uoOt/2) 

■T. [' m') dn’ (1) 

where ds is an interval along the direction of energy 
flow, and a is defined as 

a S 2/3 -t- N<r (2) 

i.e., the sum of the absorption and scattering coef- 
ficients respectively. In terms of the complex dielec- 
tric constant or of the loss tangent, tan 8, 


2/3 = [2ir(t')''VXo]e"/ ^ = [27r(€')''VXo] tan 5 (3) 

E,(z) is the directional spectral emissive power in p, 
polarization, and o>o is the scattering albedo, 

uo s Na/a (4) 

Pp,'(p, /) is the normalized scattering phase func- 
tion applicable where there exists rotational sym- 
metry. That is, Pp,/(p, p') relates the radiation in all 
azimuths of p' to that scattered into one azimuth of p. 
To paraphrase equation 1, the change in intensity 
along ds is less by the radiation absorbed and scat- 
tered and greater by the self-emission of the material 
along Bs and greater by the radiation scattered into 
the direction of Bs from all other directions. Recast 
in a more convenient form, equation 1 becomes 

p 3/p(r, p)/3t = /,(r, p) — E,/a — (wo/2) 

■ E f ' mV,(m. p') dM' . (5) 

•'-i 

where r (= —aZ), called the optical depth, is the 
optical thickness with reference to the surface. 

Equation 5 is a perfectly general integro-differen- 
tial equation describing radiative transfer within a 
scattering and absorbing medium provided that the 
scatters are randomly distributed in a matrix whose 
properties vary slowly with respect to the wavelength, 
and that scatter interaction is not dependent upon 
the path length between scatterers. In practice, the 
latter requirement, that of far-field interactions, is 
seldom realized in solid body geophysics. This con- 
stitutes a basic limitation to the quantitative, un- 
qualified application of this theory. The notation 
and formulation leading to equation 5 follow closely 
that of Chandrasekhar [1960] as will the means to 
a solution. The problem differs from most problems 
of atmospheres as treated by Chandrasekhar in that 
the dielectric is strongly self-emitting, that the scat- 
terers do not emit or absorb, and that the interface 
serves as a partial, specular reflector. One should 
note, however, that the validity of the approach is 
not compromised by these differences. 

A boundary condition for solution of equation 5 
is that the Fresnel reflection coefficients govern the 
intensity of radiation reflected into the medium by 
the surface, i.e., 

/,(0, -p) = /?p(p)/,(0, m) (p > 0) (6) 

where Rp(p) are the Fresnel coefficients. The coef- 
ficients, expressed in terms of are [Stratton, 1941) 
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iJvOx) = 


Rdix) = 


/ m - 

lAt + (t'y'w - €' + t'p.y'^1 

- (1 - + eV)'^'^ V 

+ (1 - e' + e’lyy'l 


(7) 


for (j>' less than the angle of total reflection and Rp(fj.) 
= 1 for tfi' greater than the angle of total reflection. 

The emissive power of a dielectric is dependent 
upon the complex dielectric constant and upon ther- 
mal tempjcrature, and, for isotropic media, Ls not 
dependent upon polarization. Consider an Z:(z) that 
is a linear function of z. In the absence of scatterers, 
i.e., ojo = 0, the solution to equation 5 for jj. > 0 is 


/.(r, = (Bn/2m^ ~ (E'/20E.>Xr + M>] (8) 

where 

En - E„(0) 

r - [dEM/^^l-0 (9) 

Note that equation 8 is not the general solution to 
equation 5. A particular solution to the homogeneous 
equation is Ae’’'". This, however, would apply to the 
radiation reflected downward by the surface, i,e„ 
= 0 for |n > 0 and A = /j,(0, /i) for < 0, 
Immersed in a blackbody cavity of thermal tem- 
perature Tn, the dielectric will acquire uniformly the 
temperature To, and E' will be zero. By Kirchoff’s 
law and Howell, 1968], the radiation ab- 

sorbed by the dielectric must be balanced by escap- 
ing radiation Accordingly, 

[I ~ 


The solution to the transfer equation is greatly 
simplified by assuming isotropic scattering. That is, 
radiation impinging upon a scattcrer is redistributed 
evenly in all directions. If, for further mathematical 
tractability, the energy is partitioned according to 
polarization, then 


Pp/iti, /t') = 5., (H) 

where Sp,, is the Kronecker delta. Making this sim- 
plification is troublesome because all realistic scat- 
tering phase functions have off-diagonal terms. The 
consequences of this simplification will be apparent. 
Incorporating this, the emissive power (equation 
12), and an assumed constant temperature gradient 
into equation 5 yields 

/I a/p{T, ti)/dT = /p(r, fi) — {2EdCT„/a) 

■(1 - T'r/aT,,) - (w„/2) J ' Ur, fi') d{i' (15) 

Equation 15 is not amenable to solution in'closed 
form. Physically, /pfr, /t) is expected to be continu- 
ous through T and Consequently, the integral can 
be replaced by a Gaussian quadrature, i.e., an ap- 
propriately weighted sum over 2« intervals between 
the 2 h zeroes of the even-order Legendre polynomial 
The solution achieved is at least as good as 
a power series in /j. to degree (4« — 1). The ap- 
proximation in terms of ft(= cos <}>') is appropriate 
because the intensity is symmetric with respect to 
i.e., Ip(z, +<!>') = /p(z, - «/»'). The resulting approxi- 
mation to equation 15 is 


= 111 - i?»]A'!4C0,in)...p (m > 0) (10) 

where [1 — is the fraction of blackbody inten- 
sity, hXTu), absorbed and [1 - i?„(M)] UO, 
is the intensity of thermal radiation emitted. The e' 
divisor results from the divergence as a beam passes 
into a medium of lower dielectric constant. The 
microwave blackbody intensity given by the Rayleigh- 
Jeans law [Slater, 1 939] is 

hXT)=Cr (II) 

where C is a constant. Therefore, the emissive power 
resulting from equations 8, 10, and 11 is 

E = (2j8)(e')CT (12) 

and the intensity within the dielectric becomes 

Ur> m)«.-d = e'CTotl ~ (T'/2m[r ri- (H) 

where To is the surface temperature and 7' is the 
temperature gradient. 


(x, 3/p(r, = /„(t, lii) — {20dCTo/a) 

±n 

-(1 - rr/ccT,) - (tzJo/2) E ^i,) (16) 

i-il 

where / = ±1, ■ • *, and the Oj, called the Chri'^t- 
offel numbers, arc the weighting coefficients [Love, 
1966], The zeroes and coefficients obey the follow- 
ing: 

uj = a-,;/!/ = 2 Of = 1 (17) 

J 

The solution to equation 16 is 

ft 

/p(r, Hi) — E Lpie‘‘'V(I + 

i-l 

-b «'C7-otl - T'ir + H.)/aT„] (18) 

where k) arc the n positive roots of the characteristic 
equation 
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TABLE I. Combinations of the parameters for which bright- 
ness temperatures were computed. 


wa (scattering albedo) 
e' (dielectric constant) 
a (dielectric loss) 

BT/Sz (negative thermal 
gradient) 

To (surface thermal tenapcra- 
ttire) 

tt (order of approximation) 


0.1, 0.3, 0.5, 0.7, 0.9 
1.0, 2,25, 6.25, 12,25 
1.0 (m-i) 

0, -1.0, -10.0, - lOO.OfK m'*) 
300 K 

1, 3, 5, 7, 8, 9, 10, 11. 14, 15 


1 - 0,0 E OiRy - (19) 

>>•1 

and the Lpi are defined by the boundary conditions 
at T = 0 (equation 6). Using equation 17 


/.(O. = E ^./(l 4- kill,) 

+ e^Cro[l - riii/aT^l 

4(0, —Mi) = E '^P»/(l ~ kill,) 

i-1 


+ e'CroU -h rm/aTiA (20) 



Fig, 2. Convergence of the solutions for emission dark- 
ening as a function of the order of the approximation. 
Cases are shown for an extreme scattering albedo, hjo = 
0.9, and an intermediate scattering albedo, = 0.5. 


These and equation 6 have the solution 

R„lp + - 0 <2I) 

where the elements of the matrices Rp and Ap are 

p _ / [I - ■Rn(^i)] 4- kitiA\ 4- RM]\ 

Rpi, I (22) 

and 

^ e'Cro{[l - RM)l 4- iT%/ocTo)[l + i!„(Mi)l} 

(23) 

Therefore, the elements, Lpj, of the coefficient matrix 
Lp are 


= -Efi^r'^p. (24) 

i«*L 

where are the elements of the n X n inverse 
matrix of Rp. 

The radiation intensity above the dielectric is 

4(+z, 0i) = i II - Rp(M.)lA'} MO. Ill) (25) 

where the view angle, and m, are related by Snell’s 
law, 

sin<i>i = sin (cos"‘ Mi) (26) 

7p(+z, <pi)/C has the dimension of temperature and 
is commonly called the brightness temperature, 
7’np(<#'i). Therefore, equations 18 and 25 specify the 
brightness temperature 

- {[1 - 

‘{E M.7C(1 + kiHi) + €'ro(I - T'm/aTo^ (27) 

Similarly, from equations 13 and 25, rapC^i)„,„o is 

= [1 - i?p(Mi)][ra(l - rM./2^ro)] (28) 
and the difference, caused by scattering is 

ATaMi) = TnM - (29) 

ATbp is always negative and can be called the scatter- 
induced emission darkening. 

RESULTS AND DISCUSSION 

A set of solutions for ATup (equation 29) was 
generated at each level of approximation « = 1,3, 
5, 7, 8, 9, 10, 11, 14, and 15. Each set consisted of 
solutions for every combination of parameters listed 
in Table 1. The convergence of the solutions with 
higher level approximations is illustrated in Figure 
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2. The relative convergence is more rapid the lower 
the dielectric constant. This can be understood from 
Table 2. The transmissivity of the dielectric’s sur- 
face is relatively high for ail values of m between 
p. = 1.0 and p = pcriti^l^ The transmissivity rap- 
idly drops to zero at yitpruf total reOeetton) and re- 
mains zero for a < PeritU')- The higher the dielec- 
tric constant the nearer y^ruU') is to 1,0. In the 
cases for e = 2.25, the intervals of the Gaussian 
quadrature are well-spaced with respect to pent 
(2.25). For = 12.25, however, there are only two 
intervals in the transmission window (1.0 > pi > 
pi:rit (12,25)) for n — 10, three for n — 15, and 
4 for « = 21. The fidelity of the partitioned, isotropic 
scattering model probably does not warrant going to 
the twenty-first approximation. Data in this paper 
will be from the fifteenth approximation. Fortunately, 
the more interesting AT n occur at low dielectric con- 
stants. 

Solutions to equation 29 (Figure 3) for a dielec- 
tric whose temperature is uniform (T' = 0) show 
three relationships: (a) darkening caused by scat- 
tering depends upon the dominance of scattering 
over absorption (upon (uo); (b) darkening is not ex- 
plicitly dependent upon loss, a; and (c) scatterers 
have greater effect, the lower the dielectric constant. 
Consider, for example, a homogeneous dielectric 
occupying the halfspace in Figure 1. If its loss tan- 
gent lies around 10"®, most of the emergent micro- 
wave energy will originate many wavelengths below 
the surface. Adding a few percent liquid water to the 
medium scarcely affects its dielectric constant but 
increases its loss tangent by an order of magnitude. 
Because the emissivities of such relatively low loss 
homogeneous media depend only upon their dielec- 
tric constants, adding a few percent moisture would 
not affect the brightness temperature. The source of 
the radiation, however, would tend to be far deeper 
in the dry, lower loss case. If tlie dielectric contained 
scatterers, whether they are voids or are dielectrically 
contrasting inclusions, the effect upon brightness 
temperature would be greater for the dry medium. 
In the dry case, the deeper source of the radiation 
would insure that more energy is scattered and, 
subsequently, reabsorbed. In a low loss, scattering 
medium, the introduction of small amounts of liquid 
water reduces or even effectively eliminates scaitcr- 
induced darkening. Therefore, scatterers tend to 
play an insignificant rule in emission from moist soil 
or from wet snow. In contrast, they can influence 
greatly the emission characteristics of porous ice, of 


TABLE 2. Ratio of the number of zeroes of the wlh-degrce 
Legendre polynomial within the transmission window, 1.0 > 
Mi > ficritW), to the number at angles of total reflection, in < 
lieritW)- Note that n„ii(2.25) = 0.74, ;ir„i(6.25) = 0.92, and 
;<cr.<a2.25) = 0.96. 


Order of 
approximation, n 

= 2.25 

ratio at 
e' = 6.25 

= 12.25 

1 

0/1 

0/1 

0/1 

3 

1/2 

1/2 

0/3 

5 

2/3 

1/4 

1/4 

7 

3/4 

2/5 

1/6 

8 

4/4 

2/6 

1/7 

9 

4/5 

2/7 

1/8 

10 

5/5 

2/8 

2/8 

11 

5/6 

3/8 

2/9 

13 

6/7 

3/10 

2/11 

15 

7/8 

4/11 

3/12 

17 

8/9 

4/13 

3/14 

19 

9/10 

5/14 

3/16 
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5/16 

4/17 
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Fig. 3. The effect of scattering albedo Upon emission 
darkening. is the real pari of tlic dielectric constant 
within the scattering volume. Note that darkening increases 
with the dom[n,*mce of scattering over absorption (with 
M,), and that darkening increases with decreasing dielectric 
constant. 
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Fig. 4. Emission darkening caused by scattering within a 
medium where the temperature gradient is nonzero, t' = 
2.25. Note that +j is upward and so dT/bz < 0 represents 
an increasing temperature with depth. The degree of dark- 
ening increases with decreasing dielectric constant, as is 
evident in comparisons between this figure and Figures 5 
and 6. 

dry snow, of very dry soils such as those found on 
the Moon or in some deserts, and of dry, frothy 
volcanics. 

The penetrating power of radiation in a dielectric 
depends upon the reciprocal of the loss, n. Decreas- 
ing a affects the brightness temperature by permitting 
radiation from greater depth to reach the surface. If 
the thermal temperature at depth differs from that 
at the surface, the emitted radiation will mirror that 
difference. Adding scattcrers reduces the penetration 
power. Tliercfore, an anticipated increase in bright- 
ness temperature caused by a positive temperature 
gradient with depth might well be masked by scat- 
terers. The extent of the masking is demonstrated 
in Figures 4, 5, and 6. 


Radiation emitted from the specular surface of a 
homogeneous dielectric follows the law of darken- 
ing [1 — Rp(/r)]. An integration of this directional 
emissivity over 1 > /i. > 0 yields the hemispherical 
emissivity, for each polarization, p. is slightly 
less than Sy and, therefore, the hemispherical reflec- 
tivity, (1 — ep), is slightly greater for //- than for 
I'-poIarized radiation. Radiation originating within 
the dielectric is reflected by the surface according to 
(1 — ?p) and is reabsorbed by the dielectric. Scat- 
terers, however, return some of the reflected radia- 
tion to the surface. Because reflected //-polarized 
radiation is slightly more intense than reflected 
K-polarized radiation, more H- than 1^-polarized 
radiation is returned to the surface causing AT/,/, to 
be slightly less than \Tbv- The process is apparent in 
Figure 7. 

The dependence, evident in Figure 7, of A7,,p 
upon view angle derives principally from the direc- 
tional emissivity [I — /?p(/t)]. Reducing A7«p by 



Fig. 5. Same as Figure 4, «’ = 6.25. 
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[1 ~ RpM] reveals scatter-induced changes with 
polarization and view angle (Figure 8) . Realistically, 
ATflr and should be equal at = 0. The dif- 
ference results from the assumed but artificial parti- 
tion between polarizations, Le., the scattering phase 
function, Ppii' (ji, /), should have had off-diagonal 
terms. Because the hemispherical reflectivity, (1 — 
fip), is normally greater in H- than in F-polarization, 
more scattered radiant energy is stored in the H- 
than in the F-polarized radiation fields. Without off- 
diagonal terms in the scattering phase matrix, there 
is no opportunity for a natural adjustment and, con- 
sequently, the aT/ii. differ anomalously. This anom- 
aly disappears when the hemispherical reflectivities 
are equal, as illustrated by the curve, e = 1.0, in 
Figure 8. The inference is, however, that iAT/,v| will 
tend to be greater than lAT,,,,!, and that the darken- 
ing in both polarizations will tend to be greater at 
larger view angles, 

CONCLUSIONS 

Scattering theory previously applied to atmo- 
spheric radiative transfer or to neutron scattering 



Fig. 6. Same as Figure 4, r' ~ 12.25, 



Fig. 7. Emission darkening caused by volume scattering 
as a function of view angle. Cases for two media arc shown, 
one whose dielectric constant is 2.25 and another whose 
dielectric constant is 6.25. The general shape of the dark- 
ening law, derives from the directional emtssivity 

of a homoseneoiis halfspace, [I — 7 Jj>(m)1. 

can be modified to describe the volume scattering of 
microwave radiation within dielectrics. The limita- 
tions to this particular application result from the 
assumptions of partitioned, isotropic scattering, of 
randomly distributed scattcrers, of far-field inter- 
actions (or, equivalently, of point scattering), and 
of a linear temperature variation with depth. Re- 
sulting salient inferences are; that the darkening can 
be many tens of degrees, that volume scattcrers play 
a greater role the lower the dielectric con.stant, tiiat 
the degree of darkening depends upon the dominance 
of scattering over absorption, tliat volume scattcrers 
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Fig. 8. The decrease in radiant energy available for 
emission as a function of view angle. The darkening, un- 
masked of its dependence upon directional emissivity, 
[1 — shows that relatively less energy is reach- 

ing the surface at the higher view angles. The effect is 
greater the lower the dielectric constant. Ar«i-(0) ^ 
ATirHfO) is an artifact of the partition between polariza- 
tions within the scattering phase function. The anomaly 
diminishes with decreasing dielectric constant and disap- 
pears at £* = 1.0. 


tend to mask the effects of a temperature gradient, 
that darkening tends to be greater for V- than H- 
polarized radiation, and that darkening increases 
with view angle. This solution, equation 27, rep- 
resents the functional relationships between the 
brightness temperature, the concentration of scat- 
terers, the complex dielectric constant, the surface 
temperature, and the thermal gradient. 
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DISCRIMINA.TT.ON OP GEOLOGIC MA.TERIALS 
USING SKSLAB S-192 DATA 
by Howard A. Pohn 


Introduction 

The Sliylah S-192 multispectral scanner has proved to he useful in 
the probable discrimination of several types of mate.rials in the 
southwestern Nevada desert. Although the extremes in DN Cdensity ntimber) 
in any single channel are small, a comparison within and among channels 
has shown t^t five different types of materials can be distinguished. 

Technique 

An S-192 multispectral scanner tape Caccession number 33-32882) 
for a portion of the U.S, Geological Survey Nevada test area Cfig. 1) 
was supplied to us by the Johnson Space Center. This tape was 
reformatted using the U.-S, Geological Survey DEC-1Q7Q so as to be 
compatible with our Optronics Photomation 1700 write/scan microdensitoineter . 
The images produced by the DEC-1070 were played back on film with the 
Optronics system using software written by the Remote Sensing Image 
Processing Group. 
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The images thus produced were compared with a geologic map of 
the area CStewart and Carlson, 1974) and two or three test sites were 
selected for each of five different materials. The location of the 
test sites is seen on figure 2. The materials included basalt, andesite, 
welded ash flow tuff, playa and salt marsh deposits. The DN values 
(density number, values of brightness in a 256-step gray scale) for each 
test site were extracted from the tape, and means and standard deviations 
were determined for each of the channels. The results are given in 
Table 1. 

Analysis 

It can be seen that although, in most channels, the basalts 
cannot be distinguished from andesites, channels 1 and 17 CO *52-0. 56 pm 
and 12,0-13.0 pm respectively) give a reasonably unambiguous set of 
densities which can be used to discriminate basalt from andesite. This 
observation is borne out by examining the laboratory spectral curves 
for the two materials in Hunt et al (1973, 1974). The curves show a 
relatively high reflectivity for andesite at 0,5-0. 6 pm when compared 
with basalt. This is likely due to the higher percentage of ferric iron 
and concomitant reddish color observed on the surface of many andesites. 
The ambiguity of the rest of the characteristic spectra of andesites 
and basalts may well be because of the musually low albedo of the 
andesites in this region (L. Rowan, personal communication). The higher 
DN for andesites in the daytime image at 12,0-13.0 pm is at present 
unexplained, but may be due to the lower percentage of mafic minerals 
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(and lower density) in the andesites. This would also cause the andesites 
to have a lower thermal inertia and suggests that they could be 
discriminated from basalts with thermal images. 

Welded ash flow tuffs are easily distinguishable from basalts and 
andesites from the higher DN values of the tuffs at all wavelengths 
except for 1.09-1.19 ym where it is ambiguous. 

A rather peculiar phenomenon occurs in the spectra of the playas 
and salt marshes. It can be seen that the spectra for playas are 
relatively higher than for salt marshes at all wavelengths from 
0,46-1.03 ym. However, in the Kto channels from 1.09-1.75 ym this 
relationship is reversed. At still longer wavelengths the higher DN 
values of the playas persist. 

The cause of this phenomenon is not known hut the observations 
appear to be contrary to the expected result. The channel from. 

1.55-1.75 ym occurs at the edge of a water absorption band and it would 
be expected that a salt marsh Onost probably a deposit of evaporites 
and muds containing a large amount of water) would have a lower DN 
value than a playa, which, in the Nevada desert, is typically a deposit 
of silts and clays td.th very low water content. 

Qualifications . 

The preceding statements must be qualified. First, the test areas 
have not been carefully examined in the field. Thus, it is difficult 
to say with certainty that we are detecting differences in playa and 
salt marsh deposits, because the compositions of salt marsh and playa 
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were Inferred. Second, a sampling of two or tliree areas on. a single 
SkylaL image is not felt to be sufficiently representative to make a 
broad unequivocal statement as to the capability to discriminate 
materials using spectral information. And third, the calibrated Sky lab 
S192 radiance values need to be correlated witk field and laboratory 
spectra of the test areas to see if they are comparable. 


Figure 1 Location map o£ the U.S. Geological Survey 


Nevada test area 

Figure 2. Negative print o£ S-192 channel 5 (.62-. 67 pm) 

showing location o£ the test sites. B=basalt, 
A=andesite, W=welded ash flow tuff, P=playa 
S=salt marsh. Note: the incompatibility 

between figures 1 and 2 is due to the fact that 
the S-192 tape received by the USGS was not 
corrected for the circular scan mode used on 
hoard the spacecraft. There is extreme distortion 
in the southwest; for this reason, no scale is 
provided. 

Table 1 Means and standard deviations in.DN for the test 
sites. (0 is dark, 255 is bright). 
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